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RADIOSOUNDING  OF  ICE 


V.  V.  Bogorodskiy 


Gidrometeoizdat , Leningrad,  1975 


The  physical  fundamentals  of  ice  sounding  based  on  the  pulse 
radar  method  are  presented  in  concise  form  which  is  accessible  for 
a wide  range  of  readers,  and  the  electromagnetic  characteristics 
of  various  ices  and  te^  techniques  for  their  measurements  are  also 
presented.  Also  discussed  are  the  fundamentals  of  pulse  radar  use 
for  sounding  such  comparatively  strongly  absorbing  media  as  sand 
and  permafrost,  including  permafrost  soils  and  freshwater  bodies. 
The  volume  is  intended  for  oceanologists,  hydrologists,  and  other 
specialists  involved  with  problems  of  ice  physics. 
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FOREWORD 


Man  Is  studying  the  earth  on  ever  increasing  scales.  The  con- 
siderable successes  achieved  in  this  direction  in  recent  decades  have 
led  to  the  formation  of  new  branches  of  science  — atmosphere  physics, 
ocean  physics,  earth  physics.  At  the  same  time  ice,  which  is  found 
in  large  quantities  in  the  liquid,  solid,  and  gaseous  shells  of  our 
planet,  has  essentially  been  ignored  by  the  earth  physical  sciences. 
Ice  is  studied  in  glaciology  and  cryology  — branches  of  science 
which  use  in  considerable  measure  the  geographical  and  descriptive 
methods.  However,  the  role  of  ice,  comprising  entire  systems,  is  very 
significant  in  the  geophysical  phenomena  taking  place  on  the  planet. 

At  the  present  time  the  attention  to  and  direct  interest  in  ice  -are 
constantly  increasing  in  connection  with  the  development  of  the  cir- 
cumpolar regions,  the  arctic  sea  shelf,  the  Arctic  Ocean,  and  Antarc- 
tica. Ice  has  now  become  the  object  of  broad  study  by  the  specialists 
not  only  of  those  countries  which  have  territories  adjacent  to  the 
polar  regions,  but  also  of  many  other  countries. 

The  increased  interest  in  ice,  both  as  a physical  body  and  as 
ice  sheets  (complex  geophysical  systems)  is  easily  explained.  Ice 
is  a remarkable  model  of  a proton  dielectric  — a viscoelastic  body 
which  upon  comparatively  small  changes  of  the  external  conditions 
alters  its  micro-  and  macro-structure  and  its  electrical  and  mechan- 
ical properties.  Ice  is  both  an  obstacle  for  navigation  on  the  seas, 
oceans,  rivers,  and  a structural  material.  The  continental  and  mount- 
ain glaciers  are  sources  of  pure  water  and  are  the  basic  water  re- 
sources for  irrigating  the  arid  lands  of  Central  Asia.  The  glaciers 
and  ice  sheets  are  in  considerable  degree  regulators  of  the  weather 
and  climate. 
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We  can  also  mention  the  influence  of  ice  on  life  on  our  planet 
and  practical  activity  in  man.  For  our  country,  which  as  eloquently 
expressed  by  Mendeleev,  "faces  the  north,"  the  importance  of  integrat- 
ed study  of  ice  is  difficult  to  overestimate.  In  recent  years  Soviet 
and  foreign  scientists  have  carried  out  extensive  laboratory  and 
field  studies  of  ice.  We  can  consider  that  these  studies  form  the 
basis  of  ice  physics. 

We  shall  examine  one  of  the  important  technical  advances  based 
on  the  latest  achievements  in  ice  physics  — radiosounding  of  ice. 

The  successful  resolution  of  this  complex  problem  became  possible  as 
a consequence  of  the  accumulation  of  a large  amount  of  data  on  the 
electromagnetic  characteristics  of  various  ices.  At  the  present  time 
radiosounding  of  ice  based  on  use  of  the  pulse  radar  method  has  been 
introduced  into  practice  in  the  study  of  the  antarctic  and  arctic 
glaciers  and  the  temperature  latitude  glaciers.  The  question  of  the 
radiosounding  of  drifting  sea  ice  has  been  resolved  in  principle. 

The  thickness  and  rate  of  movement  of  the  sheet-type  glaciers 
and  also  their  effective  temperature  and  the  mass  of  water  concen- 
trated in  them  can  be  measured  very  efficiently  by  the  radar  method. 

It  is  clear  from  foreign  publications  and  the  information  obtained 
during  contacts  with  scientists  of  the  U.S.A.,  Canada,  the  Scandanav- 
ian  and  other  countries  that  this  method  is  being  used  more  and  more 
widely  in  studying  the  natural  medium.  This  was  confirmed  at  the 
First  International  Symposium  on  Physical  Methods  of  Studying- Ice  and 
Snow,  held  in  Leningrad  in  October,  1973* 

The  book  presents  basically  the  results  of  studies  by  the  Arctic 
and  Antarctic  Scientific-Research  Institute,  performed  by  Trepov, 
Tripol 'nikov,  Fedorov,  Khokhlov  and  others.  We  have  also  used  the 
materials  of  certain  other  Soviet  and  foreign  studies  conducted  in 
recent  years  in  the  field  of  radioglaciology. 
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ialists involved  with  the  problems  of  ice  physics. 
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ICE  AS  AH  OBJECT  OF  GEOPHYSICAL  STUDY 


Ice,  one  of  the  aggregative  states  of  water  and  the  product  of 
metamorphism  of  snow,  is  very  widely  distributed  on  our  planet.  The 
total  volume  occupied  on  the. earth  by  glacial  ice  amounts  to  about 
25  to  30  million  cubic  kilometers,  with  nearly  90%  of  it  concentrated 
in  Antarctica.  This  is  a truly  astronomical  quantity:  there  is 

enough  ice  to  cover  the  entire  earth  with  a layer  about  60  meters 
thick.  Colossal  masses  of  ice  are  also  present  in  the  permafrost 
zone,  whose  overall  area  amounts  to  nearly  20  million  square  kilo- 
meters. The  drifting  ice  cover  of  the  Arctic  Ocean  contains  about 

60.000  cubic  kilometers  of  ice. 

Approximately  60  million  square  kilometers  of  the  land  surface 
periodically  experiences  strong  cooling,  as  a result  of  which  the 
ground  waters  lying  close  to  the  surface  and  also  the  river  and  lake 
waters  freeze,  forming  approximately  10,000  cubic  kilometers  of  ice. 

Tremendous  masses  of  pure  fresh  water  are  concentrated  in  the 

glaciers.  These  supplies  exceed  by  30  times  the  mass  of  ail  the 

water  resources  of  the  land  surface  and  by  more  than  40  times  the 

amount  of  annual  precipitation.  In  the  U.S.S.R.,  glaciers  occupy 
2 

75.000  KM  , about  one-third  of  this  area  lies  under  mountain  glaciers. 
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More  than  10  cubic  kilometers  of  water  are  present  in  the  glaciers 
of  the  U.S.S.R.  This  mass  of  water  exceeds  by  several  times  the  an- 
nual runoff  of  all  our  rivers.  The  mountain  glaciers  of  the  arid 
regions  are  very  important  for  agriculture.  Their  intense  melting  in 
the  dry  seasons  and  in  years  when  the  need  for  water  is  particularly 
acute  increases  markedly  the  flow  of  the  rivers  which  they  feed.  As 
glaciers  move,  they  alter  the  relief  significantly,  displacing  and 
then  depositing  tremendous  masses  of  moraine  material:  clays,  sands, 

and  boulders. 

The  ice  mantle  of  the  Arctic  is  more  dynamic  than  the  glaciers. 
The  ice  fields  drift  continuously  under  the  action  of  the  wind  and 
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currents.  The  nonuniformity  and  irregularity  of  the  forces  acting 
on  these  ice  fields  lead  to  the  occurrence  in  the  ice  fields  of  a 
complex  stress  state'  which  may  reach  the  critical  stress,  after  which 
the  thick  ice  fields  begin  to  crack,  break  up,  and  form  hummocks. 

The  drifting  ice  fields,  particularly  in  the  periods  of  their  com- 
pression, are  the  primary  obstacle  to  navigation  along  the  Northern 
Seaway . 

Ice  — both  as  an  element  and  as  a physical  body  — is  full  of 
contradictions.  It  conceals  an  entire  continent  from  study  and 
utilization,  interrupts  or  interferes  markedly  with  navigation, 
cements  the  soil,  destroys  coasts,  ports,  bridges.  Man  is  forced 
to  combat  ice.  On  the  other  hand,  ice  is  a source  of  fresh  water 
and  a useful  constructional  material,  from  which  nature  creates  in- 
expensive, reliable  and  comparatively  long-lived  roads,  river  cros- 
sings, takeoff  and  landing  strips,  and  floating  islands  for  drifting 
scientific-research  stations.  Ice  serves  mankind  broadly. 

We  note  that  ice  formation  and  retention  of  its  mass  involve  the 
expenditure  of  a colossal  amount  of  energy.  Thus,  the  transformation 
of  one  cubic  kilometer  of  water  at  temperature  +5°C  into  ice  at  temp- 
erature -10°C  under  ideal  thermostatic  conditions  requires  nearly  the 
same  amount  of  energy  provided  by  all  the  electric  power  stations  of 
such  developed  capitalistic  countries  as  France  or  Italy  in  the  course 
of  a year.  The  gigantic  energy  capacity  of  the  ice  formation  process 
poses  an  important  problem  of  global  nature  associated  with  study  of 
the  changes  of  the  ice  mass  on  the  planet,  since  this  ice  mass  is  both 
a product  of  the  climate  and  an  indicator  of  climatic  changes  on  the 
earth. 

2.  PHYSICAL  METHODS  OF  ESTIMATING  ICE  MASS 

The  practical  use  of  ice  in  various  economic  activities,  the 
estimation  of  the  degree  of  its  influence  on  the  geophysical  processes 
and  effect  on  various  structures,  and  the  calculation  of  the  amount  of 
water  accumulated  in  the  ice  lead  to  the  problem  of  measuring  the  ice 
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mass,  evaluating  its  state  and  rate  of  movement.  In  order  to  calcu- 
late the  mass  we  need  to  know  the  extent  of  the  ice  and  its  thickness. 
The  determination  of  the  ice  layer  thickness  obviously  causes  the  most 
difficulty . 

Several  physical  methods  have  been  used  in  recent  years  to  measure 
sheet  glacier  thickness:  gravimetric,  magnetometric,  seismoacoustic , 

and  vertical  electrosounding.  The  seismoacoustic  method  has  found 
the  broadest  application.  This  method  was  first  used  by  Matthes  in 
1926  to  measure  the  thickness  of  alpine  glaciers.  The  studies  of 
Matthes  created  the  basis  for  seismosounding  of  glaciers.  The  seis- 
moacoustic method  makes  it  possible  to  obtain  glacier  thickness  [24] 
by  direct  measurements  of  the  velocity  and  time  of  signal  propagation 
in  the  body  of  the  glacier. 

The  physical  basis  of  the  seismoacoustic  method  lies  in  reflec- 
tion of  the  sound  wave  from  the  interface  of  two  media  having  differ- 
ent acoustical  characteristics.  The  sound  waves,  usually  excited  by 
an  explosion  in  the  medium  under  study,  propagate  from  the  disturbance 
source  in  the  form  of  spheres.  As  they  encounter  another  medium 
(the  glacier  bed,  for  example),  the  spherical  waves  reflect  from  this 
medium  and  after  some  time  interval  are  recorded  by  the  receiver. 

The  time  interval  t between  radiation  and  reception  of  the  acoustic 
pulse  is  proportional  to  twice  the  distance  traveled  by  the  pulse  in 
the  medium  under  study  and  its  propagation  velocity. 


In  the  simplest  case  of  plane-parallel  layer  thickness  measure- 
ment the  sound  pulse  source  and  receiver  are  separated;  the  formula 
for  determining  the  ice  layer  thickness  will  have  the  form: 


(i) 


where  c = explosion  pulse  propagation  speed  in  the  ice;  t = time 
interval  (delay)  between  radiation  and  reception;  D = distance  be- 
tween radiator  and  receiver.  If  D is  small  [1]  simplifies  and  the 
measured  ice  thickness  will  be  determined  only  by  the  quantities  c 
and  t . 
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Seismoacoustic  glacier  studies  have  continued  for  about  50  years. 
The  seismic  method  was  first  used  in  Antarctica  in  the  fifties.  A 
major  role  was  played  by  G.  K.  Robin  in  this  work.  In  the  last  15  to 
17  years,  since  the  Norwegian,  British,  Swedish  and  French  polar 
expeditions  of  1950,  the  investigators  of  the  U.S.S.R.,  U.S.A.,  France 
Great  Britain,  Australia,  P .gium,  and  Japan  have  conducted  major 
studies  to  measure  glacier  thicknesses  and  determine  the  nature  of  the 
Antarctic  continent  rock  bed  relief. 

By  1970  seismoacoustic  sounding  of  the  glacial  mantle  had  been 
conducted  at  more  than  1000  points  along  50,000  KM-long  profiles.  A 
map  of  the  Antarctic  glacial  mantle  thickness  was  drawn  on  the  basis 
of  the  results  obtained.  The  maximal  ice  thickness  (4000  meters)  was 
noted  on  the  two  largest  Antarctic  plains:  East  and  West.  The  total 

Antarctic  ice  volume  was  established  as  24  million  cubic  kilometers, 
which  is  equivalen  to  56-meter  rise  of  the  world  ocean  level  if  the 
Antarctic  ice  were  to  melt.  The  direct  and  indirect  data  obtained 
made  it  possible  for  the  first  time  to  draw  a schematic  map  of  the 
Antarctic  sub-ice  relief.  The  mountain  massifs  discovered  beneath 
the  glacial  shield  were  termed  the  Gamburstsev,  Vernadskii,  Golitsyn, 
and  Shchukin  Mountains.  Vast  plains  were  also  discovered:  Schmidt, 

East,  and  West.  Seismoacoustic  sounding  made  it  possible  for  Soviet 
and  foreign  scientists  to  obtain  the  first  ideas  of  the  bedrock 
structure  of  Eastern  Antarctica. 

More  detailed  seismoacoustic  studies  conducted  by  foreign 
scientists  in  Greenland  also  led  to  several  scientific  discoveries. 

The  Greenland  glacier  bed  profile  and  data  on  its  thickness  are 
appearing  in  many  publications  at  the  present  time.  Soviet  special- 
ists have  made  extensive  seismic  studies  in  accord  with  the  IGY  pro- 
gram on  glaciers  of  the  Spitzbergen  Archipelago,  in  Central  Asia,  and 
other  regions. 

In  addition  to  solution  of  very  important  problems  concerning 
study  of  glacier  bed  structure  and  thickness,  the  seismoacoustic 
method  has  made  it  possible  to  gather  a large  quantity  of  data  on 
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longitudinal  and  transverse  wave  velocity  in  different  glaciers  and 
calculate  the  important  glacier  ice  physical  characteristics:  Young's 

modulus  and  shear  modulus,  and  the  Poisson  ratio.  However,  the 
primary  drawbacks  of  the  seismoacoustic  method  — impossibility  of 
remote  thickness  measurements,  low  productivity,  and  very  high  labor 
content  — have  made  it  necessary  to  seek  another  method  which  is 
free  of  these  serious  drawbacks. 

Pulse  Radar  Method  of  Studying  the  Ice  Mantle. 

In  1958  the  present  author  and  V.  N.  Rudakov  proposed  a new 
method  for  determining  glacier  thickness  based  on  the  use  of  pulse 
radar.  The  radar  (electromagnetic)  method,  just  as  the  seismoacous- 
tic method,  is  a direct  remote  method.  It  permits  measurement  of 
glacier  thickness  from  airplanes,  helicopters,  tractor-sled  trains, 
and  in  this  regard  differs  markedly  from  the  seismoacoustic  method 
and  utilizes  electromagnetic  waves  radiated  as  short  pulses  into  the 
medium  under  study  — the  glacier.  The  physical  essense  of  this 
method  is  as  follows:  the  short  electromagnetic  pulse  radiated  by 

the  antenna  from  some  carrier  vehicle  reaches  the  glacier  surface 
after  a definite  time  interval  and,  partly  reflecting  from  this 
surface,  penetrates  into  the  depth  of  the  glacier.  The  pulse  propa- 
gates in  the  body  of  the  glacier  and  after  reaching  the  bed,  reflects 
from  the  latter  and  returns  to  the  receiving  antenna.  Consequently, 
the  electromagnetic  pulse  travels  a double  distance  in  the  glacier  — 
from  the  glacier  surface  to  its  bed  and  back.  A definite  time  t is 
expended  in  traveling  this  distance.  If  we  know  the  velocity  c^  of 
electromagnetic  pulse  propagation  in  the  glacier  the  thickness  can 
be  found  from  the  formula: 


The  time  delay  t is  determined  from  the  radar  indicator  scale.  We 
note  that  the  thickness  measurement  error  does  not  exceed  2 - 2.5* 
and  the  time  of  continuous  thickness  measurement  along  unit  length 
of  the  path  being  studied  is  determined  by  the  carrier  vehicle  speed. 
Thus  the  electromagnetic  signal  delay  time  and  signal  velocity  in  the 
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body  of  the  glacier  make  it  possible  to  determine  glacier  thickness 
radpidly  and  with  high  accuracy  [5]. 

At  the  present  time,  the  radar  method  is  wisely  used  in  practice 
and  has  become  the  basic  working  "tool"  of  Soviet  and  foreign  research 
glaciologists  for  measuring  the  thickness  of  the  glaciers  of  Antarc- 
tica, Greenland,  and  other  regions.  It  is  also  being  used  with  success 
in  studying  temperature  and  mountain  glaciers  [6,  8]. 

The  pulse  radar  method,  very  effective  in  remote  measurement  of 
the  thickness  of  various  glaciers,  is  now  also  widely  used  to  study 
their  structure  and  state  (introscopy).  For  example,  this  method 
makes  it  possible  to  measure  the  average  (effective)  glacier  tempera- 
ture, determine  the  density  nonhomogeneity  and  other  nonhomogeneities 
through  its  thickness,  measure  glacier  velocity,  and  identify  crevas- 
ses hidden  by  snow  bridges.  The  radar  method  is  now  being  used  to 
measure  the  thickness  of  drifting  sea  ice  and  study  its  properties. 
Finally,  radar  glacier  studies  have  made  it  possible  to  advance 
well-founded  suggestions  on  the  possibility  of  using  the  pulse  radar 
method  for  studying  media  with  dielectric  losses  higher  than  those 
of  ice:  sands,  permafrost,  and  freshwater  bodies.  We  see  that  this 

method  makes  it  possible  to  solve  successfully  a large  range  of 
problems.  There  is  basis  to  believe  that  the  pulse  radar  method  will 
in  the  near  future  be  introduced  successfully  into  many  geophysical 
studies . 

Ice  radiosounding  is  now  based  on  extensive  theoretical  and 
experimental  studies;  they  are  discussed  in  the  following  sections. 

3.  STRUCTURE  AND  ELECTROMAGNETIC  PROPERTIES  OF  GLACIERS. 

ELECTROMAGNETIC  WAVE  PROPAGATION  IN  GLACIERS. 

Glacier  Structure. 

In  accord  with  the  Avsyuk  classification  [1]  the  largest  glaciers 
on  earth  — the  glaciers  of  Antarctica  — are  dry  polar  glaciers. 
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The  temperature  of  their  mass  is  always  below  0°C,  and  consequently 
melting  of  the  ice  ih  the  body  of  the  glaciers  does  not  take  place. 

Temperature  changes  take  place  only  in  the  surface  layer  of  thick- 
ness about  12  meters,  and  decay  below  this  level  (Figure  1). 

Glaciers  have  marked  stratified  structure.  The  surface  layer 
is  snow.  At  negative  temperatures  the  snow  layer,  primarily  as  a 
result  of  recrystallization,  transforms  into  so-called  recrystal- 

■3 

lization  firn.  The  density  of  cold  recrystallization  firn  is  0.5  g/cm 
or  higher  [25].  The  firn  crystals  have  Isometric  structure  and  their 
average  linear  dimensions  do  not  exceed  1-1.5  cm. 

The  firn  also  lies  as  a layer.  The  thickness  of  the  firn  layer 
varies,  for  example  for  Antarctic  glaciers  in  the  coastal  zone, 
the  thickness  is  about  20  meters.  In  the  cold  firn  band  the  firn 
reaches  thicknesses  of  40  meters  and  more.  Thus,  according  to  the 
data  of  the  Third  Soviet  Antarctic  Expedition  (SAE)  the  firn  layer 
thickness  in  the  Pioner  Station  region  is  120  meters,  at  the 
Komsomol ' skaya  Station  the  thickness  is  137  meters,  and  at  the  South 
Pole  the  thickness  is  164  meters.  The  firn  layer  in  turn  (for 
example,  according  to  the  data  of  Kotlyakov)  has  a stratified  struc- 
ture — more  dense  layers  alternate  with  less  dense  layers. 


The  formation  of  cold  recrystallization  ice  — the  basic  sub- 
stance of  the  Antarctic  glaciers  — is  determined  by  the  moment  of 
disappearance  of  air  permeability  of  the  firn.  It  has  been  establish- 
ed experimentally  that  the  firn  density  at  the  instance  of  air 
permeability  disappearance  approaches  0.83  g/cm  . The  glacier  ice 
layers  are  the  primary  and  thickest  glacier  layers.  Glacier  ice 
contains  certain  impurities  and  inclusions:  air  bubbles,  moraine, 

cosmic  dust,  and  so  on.  Their  influence  on  the  electromagnetic 
characteristics  of  the  ice  is  negligibly  small  in  most  cases. 

Electromagnetic  Properties  of  Ice. 

Ice,  one  of  the  basic  glacier  components,  is  a crystalline  body. 
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The  crystals  of  conventional  ice 
belong  to  the  hexagonal  group. 
The  hexagonal  ring  layer  plane 
is  termed  the  crystal  basis 
plane.  The  direction  of  the 
perpendiculars  to  the  basis 
plane  is  the  optical  axis  or 
C axis. 

Ice,  like  water,  is  a polar 
(or  dipole)  dielectric,  i.e., 
a dielectric  in  the  molecules 
of  which  the  unlike  charge 
centers  do  not  coincide  but  are 
located  at  some  distance  from 
one  another. 


In  static  fields  the  elec- 
tric properties  of  the  dielec- 
trics are  usually  described  by 


-so  -so  -<.o  -30  -n't 
o i 1 ■ 


Figure  1.  Temperature  distribu- 
tion in  upper  part  of  snow-firn 
layer  in  Antarctic  glacier  in 
’the  summer.  1-  Vostok  Station; 
2-  Komsomol ' skaya  Station;  3- 
Vostok-1  Station;  4-  Pioner 
Station. 


the  static  dielectric  permeability  ec  [3,  22].  In  alternating  fields 
the  polarization  of  the  dielectric  Is  connected  with  the  field  fre- 
quency and  also  with  the  temperature  of  the  dielectric.  In  the 
radio-frequency  band  ice  has  anomalous  dispersion,  dielectric  perme- 
ability, and  absorption.  This  phenomenon  was  explained  by  Debye, 
who  assumed  that  the  individual  molecules  rotate  under  the  action  of 
the  applied  field.  On  the  basis  of  his  theory,  Debye  obtained 
analytic  expressions,  with  the  aid  of  which  the  dielectric  properties 
of  ice  are  well  described  in  a wide  range  of  frequencies,  and  specific- 
ally in  the  ratio  frequency  band: 

‘c + «,<“->»  (3) 

1 + (<ox)2  * 


tg?— j: 

■c+‘a>(“*)2  ’ 


where  e'  and  e"  are  the  real  and  imaginary  parts  of  the  complex 
relative  to  dielectric  permeability  e (e  = e'  — je");  tg  6 — the  tangent  of 


•v 


the  dielectric  loss ‘angle;  eQ  = 93  (for  - 1 0°C)  is  the  static  dielec- 
tric permeability  of  ice;  = 3-17  is  the  dielectric  permeability 
of  ice  in  the  high  radio-frequency  band;  t = relaxation  time;  a)  = 
circular  frequency. 

The  dielectric  permeability  maintains  an  approximately  constant 
value  in  the  rf  band  and  at  low  temperatures,  while  the  dielectric 
loss  angle  tangent  decreases  with  frequency.  Both  of  these  relations 
have  been  confirmed  by  experiments,  which  makes  it  possible  to 
select  optimal  frequency  bands  and  temperatures  for  radar  system 
development.  ^ 

Electromagnetic  Properties  of  Snow  and  Flrn  in  the  RF  Band. 

Snow  and  firn  should  be  considered  a heterogeneous  two-component 
system  consisting  of  a statistical  mixture  of  dielectrics  — pure  ice 
and  air.  At  negative  temperatures  the  dielectric  parameters  of  snow 
and  firn  can  be  defined  as  the  characteristics  of  a complex  dielec- 
tric consisting  of  two  components:  pure  ice  and  air.  For  their 

determination  we  can  use  the  Lichtenecher  equations  [A]: 

= Q i®*  4 -h  ■ v 

where  • e.,,  e,. eB03a  — are,  respectively,  the' complex  relative  dielectric 
permeabilities  of  real  ice,  pure  ice,  and  air;  Q-,  and  Q2  are  the 
volumetric  concentrations  of  the  components,  where  + Q2  = 1 ; x = 
constant  quantity  for  given  distribution  of  the  components. 

Particular  cases  are  as  follows: 


1.  If  the  dielectrics  (pure  ice  and  air)  are  arranged  parallel 
to  the  vector  E,  then  x = 1 and  (5)  takes  the  form: 

• (6) 

= Qlsi.  4 “I"  Q.’-an  u* 


This  expression  must  be  used,  for  exampl'e,  when  there  are  horizon- 


tally  oriented  air  cavities  in  the  ice  and  the  electromagnetic  wave, 
polarized  horizontally,  propagates  horizontally,  or  the  vertically 
polarized  wave  propagates  vertically. 


2.  If  the  electrics  are  arranged  sequentially  relative  to  the 
vector  E,'then  x = -1  and  (5)  takes  the  form: 


l 


= Q. 


i 


~r  Q: 
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(7) 


This  will  correspond  to  horizontal  propagation  of  the  vertically 
polarized  wave  in  a glacier  with  horizontally  aligned  air  interlayers. 

3.  Finally,  the  most  typ.ical  case,  when  pure  ice  and  air  are 
"mixed"  chaotically  and  the  linear  dimensions  of  the  ice  and  air 
formations  are  considerably  less  than  the  length  of  the  electro- 
magnetic wave.  In  the  case  of  chaotic  distribution  the  component 
x - 0.  Expanding  the  terms  of  (5)  into  power-law  series,  we  obtain: 


In  *4=Qi  Ins*,. ;,-(-Q2lns*ni.  ’ 


(8) 
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For  air  in  a wide  range  of  frequencies  and  pressures  the  dielec- 
tric parameters  are  identical  to  those  for  vacuum  (e'  =1.  e"  =0). 

Since  for  ice  at  frequencies  above  10'  Hz  , Equations  (5)  - 

(8)  can  be  used  directly  for  calculating  , substituting  therein 

e'i and  . The  dependence  of  «'  on  glacier  material,  calcula- 

ted using  (8),  is  shown  for  the  213  MHz  frequency  in  Figure  2.  In 
the  calculations  we  used  the  following  basic  data:  ice  density 

p = 0.915  g/cm  , ice  relative  dielectric  permeability  e'  ,=3,355.  • 

(after  Hippel). 

Electromagnetic  Wave  Propagation  in  Glaciers. 


We  noted  above  that  the  relative  dielectric  permeability  of 
ice  is  complex: 

«-•'  ~7«W 
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The  relative  magnetic 
permeability  of  ice  is  taken 
as  one.  Considering  this  sit- 
uation and  using  the  Maxwell 
equation,  we  can  obtain  the 
fundamental  solution  of  the 
problem  of  monochromatic  wave 
propagation  in  a glacier  [18]. 


The  solution  for  the  vector 
E of  the  wave  equation  with 
complex  wavenumber  *(*)»■£.]/’«,(*), 

which  is  a function  of  the  z 
coordinate,  is  written  in 
the  form: 


(10) 


Figure  2.  r'  and  c«  as  functions 
of  density  of  glacier  material  for 
213  MHz  frequency.  1-  •i  for 
chaotic  distribution  of  components 
— pure  ice  and  air;  2-  for  chao- 
tic distribution  of  components; 

3-  for  parallel  arrangement  of 
components;  4-  for  series  ar- 
rangement; 5-  m/us;  6-  g/cm3 . 


where  A = constant. 


r 


We  consider  that  the  electromagnetic  properties  of  the  glacier 
vary  along  the  vertical  (z  axis)  because  of  the  variation  of  its 
density  and  temperature  and  write  the  complex  wavenumber  in  explicit 
form: 

*<*)««  (*)-/?(*); 1 (11) 


Here  a(z)  = propagation  coefficient: 

■ (*)— r «(*)”- r V T ( y I+i?T+ 1 ) ; 

n(z)  = refraction  index;  g(z)  = attenuation  coefficient 

•( 


pw— r*w— r 


(11a) 
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k(z)  = attenuation  index. 
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The  time  interval  t required  for  double  electronagnetic  wave 
travel  through  a glacier  of  thickness  h (before  and  after  reflection 
from  the  bed)  is  determined  from  the  expression: 


(12) 


If  the  refraction  index  varies  little  along  the  vertical  (n  = 
= const),  then 


where  ca  = electromagnetic  wave  velocity  in  the  glacier;  c = electro- 
dynamic constant. 


In  the  real  glacier  the  refraction  index  increases  with  depth, 
due  primarily  to  the  density  increase,  therefore  it  is  necessary  to 
consider  the  variation  of  n as  a function  of  the  coordinate  z.  For 
known  distribution  of  teh  density  p across  the  section  of  the  glacier 
the  integral  in  (12)  can  be  found  approximately.  The  glacier  thick- 
ness is  broken  down  into  m uniform  layers  with  refraction  index  n^ 


and  thickness  h 


k* 


r*'  * \ 

H A* 


i i 


(14) 


where  n = "effective"  refraction  index  for  the  glacier  with  known 
density  distribution  through  the  thickness. 

It  follows  from  (13)  and  (14)  that  the  time  interval  t is  deter- 
mined by  the  refraction  index  n and  thickness  of  the  glacier.  The 
effective  electromagnetic  wave  velocity  in  the  glacier: 


(15) 


Taking  for  the  frequency  f = 213  MHz  = 3-35  (after  Hip- 

pie), we  find  from  (15)  that  cr'.  = 169  m/psec. 
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Analysis  made  on  the  basis  of  (6)  and  (7)  showed  that  for  ice 
of  density  0.917  g/cm  (see  Figure  2)  the  electromagnetic  wave  velo- 
city does  not  depend  on  the  positioning  of  the  dielectrics  in  rela- 
tion to  the  direction  of  the  vector  E,  but  for  lower  densities  there 
is  some  such  dependence.  The  calculation  of  the  velocity  7*  as  a 
function  of  the  coordinate  z (and  consequently  of  p)  for  Greenland 
ice  is  shown  in  Figure  3-  It  is  clear  that  the  smaller  the  thickness 
of  the  snow-firn  layer  in  relation  to  the  overall  thickness  of  the 
glacier,  the  closer  the  quantity  c.-,  approaches  the  electromagnetic 
wave  velocity  in  pure  ice. 


The  results  of  specific  calculations  are  as  follows.  For  pure 
ice  = 3 - 35  and  cn  - 164  m/ysec.  The  thickness  of  the  snow-firn 

layer  (with  density  0.82  - 0.84  g/cm  ) for  the  Greenland  borehole  is 
73  meters.  For  the  depth  z = 200  meters  c.r  = 175  m/ysec,  for  z = 

300  meters  cn  = 171  m/ysec,  and  for  z = 400  meters  c,  - 169  m/ysec. 


It  follows  from  (10)  with  account  for  the  complex  nature  of 
the  wavenumber  that  as  the  electromagnetic  signal  propagates  in  the 
glacier  its  intensity  decreases  because  of  absorption.  The  degree 
of  signal  attenuation  (in  dB)  can  be  found  from  the  formula: 

JV,=20lg-f^  (l6) 

where  Eq  = radiated  signal  intensity;  E = received  signal  intensity. 

For  plane  waves  which  have  twice  traversed  a glacier  of  thick- 
ness h we  obtain: 

(17) 


= 1 7.36  f ? (z)  rfz—  1 7.36  f A (*)  <&. 
I e i i 


For  high  frequencies  tgrt<|.  . Assuming  constancy  of  tg  6 through 
the  thickness  of  the  glacier,  we  obtain  the  final  expression  for  N„ 

AfB-8.68~/i7/Hg8. 


n 


(18) 
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Thus,  from  theoretical  exam- 
ination of  electromagnetic  wave 
propagation  in  a glacier  we  have 
obtained  on  the  basis  of  solu- 
tion of  the  wave  equation  (10) 
the  time  lag  between  the  radia- 
ted (sounding)  and  received 
radar  signals  and  the  magnitude 
of  the  absorption.  These  quant- 
ities make  possible  quite  ac- 
curate determination  of  the 
parameters  of  radars  for  glacier 
thickness  measurement.  Electro- 
magnetic wave  propagation  in 
glaciers  is  also  accompanied  by 
other  losses.  They  will  be 
examined  later. 


Figure  3.  Distribution  of  glacier 
material  density  along  the  vertical 
and  dependence  of  on  distance 
to  the  surface  of  the  Greenland 
glacier.  1-  p g/cm3 ; 2-  m/ys. 


4.  EXPERIMENTS  TO  DETERMINE  ELECTROMAGNETIC  CHARACTERISTICS  OF 
GLACIERS 
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Measurements  of  Electromagnetic  Wave  Propagation  Velocity. 


It  follows  from  (13)  that  knowledge  of  the  electromagnetic  wave 
velocity  for  the  problem  under  consideration  is  very  important, 
since  it  determines  the  accuracy  of  the  method.  The  most  reliable 
method  for  measuring  electromagnetic  wave  velocity  is  the  slant 
sounding  method.  This  method  was  used  successfully  by  Gudmandsen 
in  Greenland  and  by  Trepov  and  Jiracek  in  Antarctica. 

The  constancy  of  the  electromagnetic  wave  propagation  velocity 
in  ice  in  a broad  frequency  range  (from  10  to  2700  MHz),  its  sm?ll 
variation  as  a function  of  temperatures  existing  in  natural  glaciers 
and  the  large  amount  of  data  on  the  structure  of  these  glaciers  make 
it  possible  to  obtain  reliable  data  on  the  value  of  c.,  for  various 
geographic  coordinates  and  experimental  locations. 
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The  essence  of  the  slant 
radio  sounding  method  is  illus- 
trated in  Figure  4.  The  trans- 
mitting antenna  with  the  radio 
pulse  transmitter  is  located  at 
point  B.  -The  receiving  antenna 
and  the  receiver  with  indicator 
are  located  at  point  A,  at  the 
distance  1_  from  the  transmitter. 
To  each  pulse  radiated  by  the 
transmitter  correspond  two  pulses 
in  the  receiver:  one  travels 


Soidyx 


A L .8 


Figure  4.  Scheme  for  determining 
electromagnetic  wave  velocity  by 
slant  sounding.  1-  air;  2-  ice; 
3-  glacier  bed. 


the  distance  1 through  the  air, 

the  other  is  reflected  from  the  glacier  bed  and  travels  the  distance 
s through  the  ice  (in  the  figure,  point  A'  is  the  mirror  image  of 
point  A relative  to  the  glacier  bed).  Measuring  on  the  Indicator 
the  delay  of  the  second  signal  relative  to  the  first  for  different  1, 
we  can  from  geometric  relations  obtain  the  velocity  of  the  radio  pulse 
in  the  ice  under  the  condition  that  the  glacier  is  homogeneous: 


ity  in  the  ice  reduces  to  obtaining  the  experimental  relation 

approximating  the  experimental  relation  by  a straight  line,  and 

determining  the  slope  of  this  approximating  line.  The  time  t is 

found  as  the  sum  of  the  delay  time  t^  of  the  pulse  traveling  through 

the  ice  relative  to  the  pulse  propagating  in  the  air  (measured  on  the 

indicator)  and  the  time  of  pulse  propagation  through  the 

c 

air.  The  distance  1 is  measured. 

The  experiments  were  conducted  as  follows.  The  mobile  receiving 
and  transmitting  sets  were  moved  in  opposite  directions  along  a pre- 
determined rectilinear  survey  line  so  that  the  reflecting  bed  segment 
remained  the  same  in  the  process  of  the  measurements.  Measurements 
were  made  every  50  meters  of  the  lag  of  the  pulse  reflected  from  the 
bed  relative  to  the  pulse  traveling  through  the  air.  Preliminary 
measurements  of  the  glacier  thickness  were  made  before  starting  the 
study  in  order  to  select  a segment  with  constant  thickness.  The 
results  of  velocity  measurements  using  a radar  with  carrier  frequency 
440  MHz,  obtained  by  the  14th  SAE  [20]  along  one  survey  line  near 
the  Molodezh  Station  are  shown  in  Table  1. 


TABLE  1 
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It  is  well  known  [16]  that  the  slope  of  the  straight  line  approx- 
imating the  experimental  linear  relation  y = a + bx  between  two  random 
quantities  x and  y when  using  the  least  squares  method  is  defined  by 
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where  r = coefficient  of  correlation  between  x and  y. 
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The  results  of  measurements  of  along  seven  survey  lines 

are  shown  in  Figure  5-  The  measured  radio  wave  velocities  167-8, 
168.6,  170  and  178  m/ysec,  obtained  near  the  Antarctic  shoreline, 
and  184-5  m/ysec  obtained  45  km  from  the  coast,  agree  well  with 
the  data  calculated  using  the  average  glacier  density.  The  experi- 
mental data  basically  confirm  the  validity  of  the  velocity  c„  = 

= 169  m/ysec  used  for  the  calculations.  Analysis  of  the  measurement 
error  showed  that  the  basic  error  apparently  occurs  in  reading  out 
the  time  interval  between  the  pulses,  since  for  a baseline  length 
approximately  equal  to  the  glacier  thickness  this  interval  is  short 
(in  our  case,  0.55  ysec)  and  is  comparable  with  its  readout  accuracy. 
Ir.  order  to  reduce  this  error  we  must  make  measurements  on  a baseline 
with  />  2 - 3h.  Such  measurements  are  possible  with  small  glacier 
thicknesses . 

In  the  case  of  long  baselines,  when  the  wave  incidence  angles 
on  the  glacier  surface  are  large,  error  may  arise  due  to  refraction 
at  the  air-glacier  interface  and  in  the  upper  layers  of  the  glacier, 
which  have  variable  coefficients  of  refraction.  In  this  case  as  a 
result  of  increase  of  the  path  length  in  the  layer  with  smaller 
refraction  coefficient  the  radio  wave  velocity  will  be  higher  than 
the  velocity  obtained  at  smaller  incidence  angles.  However,  this 
error  is  systematic  and  can  be  taken  into  account  theoretically  on 
the  basis  of  data  on  the  antenna  height  above  the  glacier  surface  and 
the  density  vertical  distribution. 

Table  2 shows  the  velocities  obtained  by  the  15th  SAE.  The 
measurements  were  made  at  a 440  MHz  frequency  at  three  points  in  the 
Molodezh  Station  region.  At  each  point  the  measurements  were  made  on 
three  different  survey  lines  with  the  use  of  parallel  and  orthogonal 
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receiving  dipoles. 

In  the  I6th  SAE  velocity 
measurements  were  made  for  the 
first  time  at  60  MHz  frequency. 
The  use  of  a powerful  radar 
set  made  it  possible  to  conduct 
measurements  in  regions  with 
large  glacier  thickness. 
Measurements  were  made  at  six 
points  on  a survey  line  to  the 
south  of  Molodezh  Station. 

The  measurement  results  are 
shown  in  Table  3. 

Jiracek  (USA),  using  a 
similar  method,  obtained  at 
the  30  MHz  frequency  the  follow- 
ing radio  wave  propagation 


Figure  5.  Experimental  data  and 
approximating  curves  of  t§  versus 
1 * . Average  values  of  1;  2- 
168.6  m/ysec;  3-  178.0  m/ysec; 

4 & 5-  184.5  m/ysec;  6-  172 . 0 m/tisec ; 
7-  167.0  m/ysec;  8-  ys. 


velocities  in  ice:  with  glacier  thickness  from  170  to  320  meters 

c.i  * 180  + 5 m/ysec,  on  the  760  m-thick  Skelton  Glacier  (without 
snow-firn  cover),  c = 168.5  m/ysec,  on  the  Roosevelt  Glacier  with 
thicknesses  500  and  800  meters  r,  = 174.8  m/ysec.  These  velocities 
make  it  possible  to  obtain  the  average  dielectric  permeability  and 
refraction  index  of  radio  waves  for  the  glaciers  located  in  the 
region  of  the  experiment. 


TA3LE  2 


Measurement  line 
identification  No. 

1 

2 

3 

1 

2 

3 

1 

2 

3 

Average  glacier 
thickness  along 
line,  m 

300 

30i) 

309 

420 

,V>8 

4M 

3fi7 

384 

375 

Radio  wave  velo- 

city  in  ice,  m/usec. 

170 

ir.r> 

I8!i 

170 

183 

179 

172 

176 

' i 22 

I 
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TABLE  3 

RESULTS  OF  RADIO  WAVE  PROPAGATION  VELOCITY  MEASUREMENT  IN  ICE 
ON  SURVEY  -LINE  TO  THE  SOUTH  OF  MOLODEZH  STATION 


Point 

Number 

Direction 

Radio  Wave 

Velocity, 

m/us 

Glacier 
Thickness , 
m 

43 

Northward 

172 

Southward 

160.5 

820 

46 

— 

167.0 

597 

22 

Northward 

182.0 

860 

Southward 

172.5 

109 

Northward 

162.5 

650 

Southward 

169 

26 

Northward 

153 

1330 

Southward 

154 

73 

Northward 

160 

313 

Southward 

170 

Study 

of  Dielectric  Losses 

in  Antarctic  Glaciers. 

It  was  mentioned  above  that  ice  is  a dielectric  having  complex 
dielectric  permeability.  The  Imaginary  part  of  the  complex  dielectric 
permeability  determines  the  magnitude  of  the  dielectric  losses.  To 
characterize  the  ability  of  a dielectric  to  absorb  electromagnetic 
field  energy,  we  generally  use  the  dielectric  loss  angle  tangent  tgo, 
which  is  defined  as  the  ratio  of  the  two  components  of  the  complex 
dielectric  permeability  or  as  the  ratio  of  the  active  conductivity 
of  the  imperfect  dielectric  to  its  capacitive  conductivity,  i.e., 

mf*  Vi 

tgd— — . The  tangent  of  6 is  a very  inraortant 

e . we  „• 

electromagnetic  characteristic  and  is  necessary  for  calculating 
electromagnetic  wave  absorption  in  ice  and  selecting  the  radar  poten- 
tial. Field  measurements  of  tg  6 are  most  interesting  for  the  subject 
problem. 


In  1966  the  range  of  tg  6 values  for  actual  Antarctic  ice  was 
calculated  indirectly  [10].  The  data  of  radar  sounding  along  the 
route  from  Mirnyi  Observatory  to  Pioner  Station  were  used  for  this 
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TABLE  4 


TANGENT  OF  DIELECTRIC  LOSS  ANGLE  OF  REAL  ICE  (AFTER  JIRACEK ) . 


Ice 

Ice 

Study  Region 

Density, 

g/cm3 

Temper- 
ature, °C 

tg5 

for  frequency,  MHz 

■ ISO 

• • 

300 

Iceland 
(glacier  ice) 

0,818 

1 • 

• Y 

•» 

-1 
i —5 

-io  . 
-20 
-ao 

—40  • 

■•-50  ■- 

—50- 

* • r *’ 

0,0022 
0,00144  '• 
0.00110 
0.00068 
0.00043 
0.00020  ‘ 
• 0.00014 
0.00006 

• i t-»  *r**»»»*r  *• 

0.00108 
0.00076  . 

o.ooors 

0.00033 

o.ooo3r 

0.00013  i 
0,00008  < 
0,00003., 

r>~!T 

0.00052  d 
^ 0.00040  — 
a 0.00088. { 
0.00019. 
0,00013  * 
o.oooood 
0,00006^ 
0,0000*,;? 

Little  America, 
Antarctica 

0.881 

* ♦ ‘ • 

• * '■/» 
1 

-1  ■' 
—5' 
10«., 
-20 
—30 
■’  —«0 

"0.0037 
•0.0086 
.0,00217  i 
0,00154 
. 0,00116  . 
0.00085 

o.ooifr 
:•  0.0013  1 
? 0.00106** 
0,00078  • 
0,00067,' 
0,00044  * 

■'ro,ooioo"J« 
>•0.00072$ 
*1*0.00068^1 
* 0.00033  d 
< 0,0008*? 
0.00085  ^ 

purpose.  After  reducing  the  measurement  data  a calculation  was  made 
of  the  dielectric  losses  resulting  from  absorption  for  various  values 
of  tg  6,  and  then  the  overall  losses  were  evaluated  for  each  assumed 
value  of  tg  6.  It  was  concluded  in  comparing  the  experimental  and 
calculated  losses  that  the  most  reliable  average  values  of  tg  6 at 

the  213  MHz  frequency  for  the  entire  thickness  of  the  glacier,  in- 

— 4 —4 

eluding  the  snow  layer,  lie  in  the  range  4.2  • 10  - 7 * 10 

Jiracek  presented  the  most  complete  data  on  the  values  of  tg  5 
for  real  ice  in  1967.  (Table  4).  We  see  from  the  table  that  tg  6 
depends  strongly  on  ice  temperature  and  field  frequency. 
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5. 


ATTENUATION  OF  RADIO  WAVES  IN  GLACIERS 


Overall  Attenuation  of  Radio  Waves  During  Their  Propagation  In 
Glaciers . 

It  was  shown  above  that  radio  pulses  decay  as  a result  of  ab- 
sorption during  passage  through  a layer  of  ice.  It  follows  from 
( 1 8 ) that  radio  wave  absorption  is  determined  by  the  radiated  fre- 
quency (the  working  frequency  f is  connected  with  the  circular  fre- 
quency cu  by  the  relation  1),-  ice  layer  thickness,  and  the  quantities 
e*  and  tg  6.  Absorption  is  an  irreversible  loss  of  radio  wave  energy; 
the  field  energy  is  transformed  into  heat  as  a result  of  interaction 
with  the  medium. 


Radio  waves  also  decay  as  a result  of  spreading  of  their  front 
during  propagation  in  the  medium  (geometric  losses)  and  reflection  on 
the  interfaces  (surface,  interior  density  nonhomogeneities,  glacier 
bed).  They  also  change  their  intensity  as  a result  of  the  focusing 
action  of  the  ice  layer  and  noncoincidence  of  the  reflected  signal 
polarization  with  the  receiving  antenna  polarization.  Study  of  the 
components  of  radar  signal  attenuation  during  propagation  in  the 
various  earth  covers  is  important  primarily  from  the  viewpoint  of 
optimizing  the  parameters  of  radars  intended  for  sounding,  and  also 
in  refining  the  electrodynamic  models  of  these  covers  and  studying 
their  electrophysical  characteristics. 


Radio  wave  propagation  theory  and  the  large  amount  of  experiment- 
al data  obtained  in  radar  sounding  of  the  Antarctic  and  Arctic 
glaciers  make  it  possible  to  evaluate  quantitatively  the  overall 
attenuation  on  the  basis  of  the  previously  indicated  additive  compon- 
ents of  electromagnetic  wave  attenuation  in  ice.  The  amplitudes  of 
the  signals  received  by  the  station  after  their  passage  through  the 
ice  layer  are  determined  by  the  overall  attenuation,  which  can  be 
expressed  by  the  relation: 


, V:  - /Vr  -f  /Vntp  • i-  ,V.„ + iVp  -r  A/„  -f  /V,„ 
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where  AV  = geometric  losses;  V„r|,  = losses  owing  to  reflecting  from 

the  interfaces;  A'*  = signal  change  as  a result  of  focusing;  A',.  = los- 
ses as  a result  of  scattering;  N„  * losses  resulting  from  absorption; 
A'm.i  = losses  resulting  from  noncoincidence  of  the  polarization  of 
the  signal  received  by  the  antenna  with  the  receiving  antenna  polar- 
ization. 


We  shall  consider  that  the  signal  reflected  by  the  glacier  dif- 
fers from  the  radiated  signal  only  in  amplitude,  time  interval,  and 
phase  shift.  This  signal  acts  together  with  noise  on  the  receiver 
input.  The  noise  level  is  determined  by  the  noise  level  of  the 
receiving  channel,  the  atmosphere,  and  the  medium  being  sounded. 

If  we  examine  the  problem  of  measuring  glacier  thickness  we  can 
consider  that  part  of  the  input  is  the  useful  signal  (in  our  case  the 
pulses  reflected  by  the  glacier  bed),  while  the  remainder,  consisting 
of  noise  and  the  pulses  reflected  by  the  upper  boundary  and  the 
internal  nonhomogeneities  of  the  glacier,  is  interference.  We  will 
not  examine  the  theory  of  received  signal  detection  and  processing; 
this  theory  is  examined  in  detail  in  many  specialized  monographs 
and  texts.  A more  important  problem  is  the  analysis  and  evaluation 
of  the  terms  determining  the  magnitude  of  the  overall  radio  signal 
attenuation  during  sounding  of  media. 


Absorption 


It  was  shown  previously  that  the  dielectric  loss  angle  tangent 
depends  strongly  on  the  temperature  of  the  medium  in  question.  If 
the  vertical  temperature  distribution  through  the  glacier  is  known  the 
absorption  can  be  calculated  using  (17),  which  takes^.nto  account  sig- 
nal propagation  toward  the  glacier  bed  and  back.  If  k(z)  is  repre- 
sented by  a complex  function  of  z the  losses  due  to  absorption  can 
be  evaluated  approximately  by  considering  the  glacier  to  consist  of 
a large  number  of  layers  with  dielectric  loss  angle  tangent  which 
is  constant  in  the  limits  of  each  layer.  The  losses  in  the  layer  are 
calculated  from  an  expression  similar  to  ( 1 8 ) : 
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N^z)- 8.68 tg  Idz 


(23) 


where  dz  denotes  a very  small  thickness  of  the  layer  with  constant 


tg  6;  Nn(z)  is  the  absorption  in  the  layer. 


Now  we  must  sum  the  losses  in  the  layers  through  the  entire 
thickness  of  the  glacier 

» (24) 


Na^N„(z)dzL 


or,  if  the  number  of  layers  is  finite  (for  example  m layers) 

m 

2/V0(*)  Az, 


(25) 


where  Az  = thickness  of  layer  with  practically  constant  tg  6.  It  is 
obvious  that  the  thinner  the  layers,  the  more  accurately  the  approx- 
imate expression  (25)  evaluates  the  losses,  and  the  closer  it  ap- 
proaches the  more  exact  expression  (24). 


However,  in  most  cases  the  temperature  distribution  through  the 
section  of  the  glacier  is  not  known.  Therefore  it  is  convenient  to 
characterize  the  absorbing  properties  of  the  glacier  body  by  the  so- 


called  absorption  temperature  t , satisfying  the  expression: 


tf.-AUO*. 


(26) 


where  Nn(t„)  is  the  specific  absorption,  uniquely  defined  by  the  ab- 
sorption temperature. 


The  absorption  temperature  thus  characterizes  the  glacier  as  a 
whole,  specifically  its  ability  to  absorb  electromagnetic  signals, 
and  is  a very  important  additional  thermal  characteristic  of  the  gla- 
cier, since  it  makes  it  possible  to  differentiate  warmer  and  cooler 
glaciers  without  reference  to  the  thickness. 


The  values  of  ,V„(/11)  were  calculated  on  the  basis  of  ( 1 8 ) for 
unit  glacier  thickness  (per  meter)  and  are  shown  in  Table  5 for  the 
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10  - 500  MHz  frequency  range. 


TABLE  5 


fn  *0 

A’„  . . . 

/ 

* II  -•••** 

A',,  . . . 


-1 

— *1 

-10 

-20 

U.07I  ft 

o.lho 

0.0122 

. . --.Vi 

— 10 

— bO 

. 0.1120  » 

0.01“*) 

0,00*7 

0.0025 

1-  dB/m 

These  data  make  it  possible  to  calculate  easily  the  radio  wave 
energy  losses  by  absorption  for  any  ice  temperature. 

Geometric  Losses  and  Losses  at  the  Lower  Boundary. 

In  accordance  with  the  radar  range  expression,  we  can  write  for 
the  geometric  losses  and  the  reflection  losses  AV  at  the  lower 
boundary  for  the  cases  of  mirror  reflection  and  scattering  with  the 
radar  set  operating  in  the  near  zone*:  ( h <20  + 30  --  ---  J 

iV,  - j- <Votp-  - ‘JO  In  <//•  t 20  In  * - | 10  Ig  64f  — 20  Ig  Kh ' 

where  GA  is  determined  by  the  station  parameters;  G = radar  antenna 
gain;  Kf  = Fresnel  reflection  coefficient;  for  media  with,  small  losses 
and  normal  incidence:  ....  . _ 

K,~  j/r1  + jA; ' 

\„rt,  =»  — 20  Ig  Kf  and  depends  on  the  dielectric  permeability  of  the 
underlying  surface.  For  the  calculation  we  took  jVotp  = 14  dB. 


Here  r.„,  = duration  of  rf  pulse  radiated  by  the  radar 


r 


Summing  the  geometric  losses  and  the  reflection  losses,  we  obtain 
( in  dB > : . AV4-  A'„r|.=~  —20  Ik  Gl-J-SO  Ig  h }-42. 

(27) 


Focusing. 

Jumping  ahead,  we  shall  show  that  signal  attenuation  as  a con- 
sequence of  focusing  when  working  from  the  glacier  surface  N » =■  — 20  I g n j 
where  n = coefficient  of  refraction  at  the  air-ice  boundary.  When 
working  from  an  airplane  it  is  necessary  to  consider  the  dependence 
of  focusing  on  the  flight  altitude. 

We  shall  examine  briefly  some  theoretical  aspects  of  electro- 
magnetic wave  focusing  when  sounding  glaciers. 


It  is  well  known  that  in  a homogeneous  medium  at  the  distance  r 
from  a point  radiation  source  the  power  flux  density: 

a,.  (28) 


//«- 


mi 

Ur- 


In  a stratified  nonhomogeneous  medium  the  power  flux  density  at 
the  same  distance  will  be  different  because  of  refraction  and  equal 
to  n.  We  term  the  ratio: 


n_ 

no 


(29) 


the  focusing  factor.  We  shall  determine  the  focusing  factor  for  a 
stratified  nonhomogeneous  medium  in  which  the  refraction  coefficient 
n is  a function  of  the  z coordinate.  The  case  of  two  media  with  a 
flat  interface  will  obviously  be  a particular  case.  In  Figure  6 the 
radiation  source  is  located  at  the  coordinate  origin,  aQ  = ray  glance 
angle  at  the  exit  point;  = ray  glance  angle  at  the  observation 
point  at  the  depth  h.  We  obtain  ds'=  2nR  sin  a*  dR:  for  the  wave  front 

area  included  between  the  rays  leaving  at  the  angles  on  and  ao  + Wa. 

dR 

Considering  that  — — da « (the  minus  sign  means  that  for  h = 



const  smaller  R correspond  to  larger  otQ),  we  obtain: 
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Figure  6.  Illustration  for  determining  focusing 
factor  for  stratified  nonhomogeneous  medium. 


’7X3 

Dn 


Since  the  power 
is  equal  to 


in  the  solid  angle  ft  in  the  limits  from  ct*'  to 
PQmm— ^^cosaoda^  , the  power  flux  density 


at  the  distance  r from  the  source  is: 

1 'rcr-r-.rp 


\ljmm  -V 


Substituting  (28)  and  (30)  into  (29),  we  obtain: 

r2  cos  io 


(30) 


(31) 


t=- 


_ d/?  , 

* "3^ sln  *» 


In  order  to  calculate  the  focusing  factor  using  (31)  it  is 

OUtl  (Iff 

necessary  to  determine  R and/?*  The  final  expression  for  n 


has  the  form: 


cos2 


‘ d£ 

10  J |«2  (.?)  — cos2  »ol ‘/* 


sin  ih  sin  iq 


n2  ( z ) ,/* 


{ |n2  (*)  — cos*’  ,o| 


(32) 
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In  the  case  of  a single  refracting  boundary  n(z)  = nQ  for 
0<*4;Ai  and  n(z)  = n^  for  h.  , where  hQ  = antenna 

above  the  glacier  surface.  Then  with  account  for  the  relation 
a.  = 1 0.,: (32)  takes  the  form: 

n ii 


hg  _4_  ^ 

It; 

h q sin  3|  h slu  io 

n0  sin- 10  n|Sln-*i 


height 

sin 

(33) 
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Table  6 shows  the  results  of  calculation  of  the  focusing  factor 

n for  radiation  passage  from  a medium  with  refractive  index  n^  = 1 

(air)  into  a medium  with  refractive  index  n,  = 1.5  (snow  with  density 

3 

0.5  g/cm  ) and  back  for  various  glance  angles  in  the  first  medium 
(air). 
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i 

\ 


\ 


! 


i 

i 

1 
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TABLE  6 


Glance  angle  in  first 
medium,  degree  . . . 


'for  n =1  and 
o 

n1  = 1.5  

f°r  = 1.5  and 

n,  = 1 

■»  i 

Attenuation  for  radiation 
propagation  in  both  direc- 
tions, dB 


1 

2 

5 

0.020 

0.0T.2 

0.13 

VI® 

- ? 0ytt.tr 



- - . . 

..  

_xj ■"s.aiiil 

0,011 

0,023 

0.055 

35,3  29.3  21.3  J 


Let  us  evaluate  focusing  of  electromagnetic  energy  at  small 
incidence  angles  (vertical  radiosounding).  For  this  case  (31)  for 
incidence  angles  fi=90°  — a • has  the  form: 


rJMn  ?n 


OR 

*V) 


CO*^ 


(3^) 


Omitting  the  intermediate  arguments,  but  considering  the  rela- 
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tions,  valid  for  small  incidence  angles. 


( os  ?o  cos  ss  \ . 
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we  obtain  the  formula  for  n in  the  plane  interface  case: 
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(35) 


The  focusing  factor  n * for  sounding  of  the  medium  is  equal  to 
the  product  of  the  focusing  factors  for  wave  propagation  "to" 
and  "from"  (norp): 


I 4 


_h 
ho  I 


V 


(36) 


where  n— is  the  relative  refraction  index  at  the  boundary  of  the 

no 

two  media. 

We  note  that  (36)  was  obtained  for  the  field  of  a point  source 
located  at  the  point  where  the  antenna  is  located  (direct  signal 
focusing)  or  at  the  point  where  the  fictitious  source  is  located  for 
the  reflection  case  (reflected  signal  focusing).  The  radiator  or 
brightness  is  assumed  to  be  the  same  in  nonhomogeneous  and  homogeneous 
media.  The  focusing  factor  n ' does  not  consider  the  variation  of 
the  effective  reflecting  surface  of  the  object  which  may  occur  with 
replacement  of  the  real  nonhomogeneous  medium  by  a homogeneous  (ref- 
erence) medium. 

We  note  also  that  the  magnitude  of  the  focusing  factor  will  de- 
pend on  the  geometric  characteristics  of  the  reflecting  surface. 

The  variation  of  the  effective  reflecting  surface  of  some  objects 
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located  in  a nonhomogeneous  medium  was  analyzed  by  Trepov.  Figure  7 

shows  curves  of  n ' as  a function  of  the  ratio  h/h  for  reflection 

o 

from  a point  target  in  accordance  with  (36),  from  a mirror  surface, 
and  from  a scattering  surface  located  in  the  far  zone  in  relation 
to  the  radar. 

Analyzing  these  relations,  we  can  say  that  in  the  case  of 
vertical  sounding  of  a medium  from  air  the  contribution  of  the  focus- 
ing factor  to  the  overall  radar  signal  attenuation  is  not  always 
significant.  In  all  cases  it  lies  in  the  limits  from  -20  lg  n to 
+ 20  lg  n,  which  for  sounding  of  ice  from  air  varies  from  -5  to 
+5  dB.  However,  there  is  no  doubt  that  the  dependence  of  the  focus- 
ing factor  on  the  ratio  h/hQ  is  one  of  the  reasons  for  the  strong 
dependence  observed  in  practice  of  radar  signal  attenuation  on  air- 
plane flight  altitude  in  radar  measurements  of  glacier  thickness. 

It  is  precisely  for  this  reason  that  glacier  sounding  is  usually 
conducted  from  low  altitudes. 


Scattering. 


Let  us  examine  radar  signal  attenuation  as  a result  of  scattering 
of  its  energy  at  the  upper  boundary  of  and  in  the  medium  being 
studied.  In  practice  two  problems  arise:  one  is  accounting  for  the 

influence  of  scattering  on  the  overall  radio  signal  attenuation  during 
sounding  of  the  medium,  the  other  is  evaluation  of  the  magnitude  of 
the  signal  "back"  scattered  by  the  surface  irregularities  and  the 
nonhomogeneities  within  the  medium,  since  this  signal  is  noise  when 
measuring  the  glacier  thickness  and  may  mask  the  useful  signal  re- 
flected from  the  bed. 


Usually  the  scattering  capability  of  the  medium  is 
by  the  scattering  index  c and  the  scattering  indicatrix 
are  defined  by  the  relations: 


df\  . / \ 
-jn1- ' *(T)= 

n m 


r(  t) 
C 


characterized 
x(y) , which 


where  p„,a  and  are  the  powers  incident  in  the  volume  dV  and 
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Figure  7.  Variation  of  focusing  factor  n ' with 
reflection  of  radio  signal  from  point  target  (a)  , 
from  mirror  surface  (b),  and  from  scattering 
surface  in  the  far  zone  (c).  1-  a,  2-  b,  3-  c. 


scattered  by  this  volume  in  all  directions;  p(y)  is  the  power  scat- 
tered in  the  unit  solid  angle  in  the  direction  forming  the  angle  y 
with  the  incident  wave  direction. 


Given  these  characteristics  and  the  absorption  coefficient  of 
the  medium,  the  "forward  and  back"  scattered  field  can  be  determined 
by  solving  the  transport  equation  [19].  These  characteristics  have 
not  yet  been  studied  in  the  rf  band  for  snow,  firn,  and  ice,  and 
therefore  exact  solutions  of  the  problems  indicated  above  are  not 
possible  at  the  present  time. 


t 
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An  approximate  estimate  of  the  attenuation  caused  by  scattering 
in  a medium  can  be  based  on  the  fact  that  the  considered  media  have 
a stratified  structure  with  approximately  parallel  layer  boundaries. 

In  certain  cases  (water  layer  in  a glacier,  thawed  soil  layer  in 
permafrost)  the  attenuation  caused  by  scattering  may  be  significant. 

The  signal  attenuation  resulting  from  scattering  at  the  upper  bound- 
ary can  be  taken  into  account  approximately  by  the  Fresnel  coefficient, 
since  for  small  incidence  angles  the  reflection  coefficient  is 
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independent  of  surface  roughness.  It  follows  from  the  general 
principles  of  wave  scattering  theory  that  the  interface  "forward" 
scattering  diagram  is  brooder  for  radiation  passage  from  the  medium 
into  the  air  than  in  the  reverse  direction. 

We  estimate  the  losses  (in  dB)  due  to  scattering  as  the  losses 
in  passage  through  the  upper  surface  of  the  glacier: 

A^-20  lg(l -K\  )«0.7, 

where  for  the  air-glacier  interface  K//"0,284.', 

It  is  obvious  that  the  losses  at  the  low-contrast  boundaries  of 
the  density  nonhomogeneities  in  the  glacier  will  be  still  smaller 
and  can  be  excluded  from  consideration. 

Polarisation  Losses. 

It  is  well  known  that  if  we  study  an  electromagnetic  wave  with 
linear  polarization  in  an  isotropic  medium  (air,  for  example),  the 
signal  received  from  a reflecting  object  will  be  maximal  when  the 
receiving  dipole  is  parallel  to  the  radiating  dipole  (provided  the 
depolarizing  properties  of  the  interfaces  are  not  significant).  In 
other  words,  an  isotropic  nongyrotropic  medium  does  not  have  de- 
polarizing properties. 

It  was  discovered  by  the  14th  SAE  that  the  signals  reflected 
from  the  glacier  bed  have  polarization  differing  from  that  of  the 
radiated  signal.  Thus,  with  linear  polarization  of  the  radiated 
signal  the  reflected  signals  have  elliptic  polarization  (linear 
in  particular  cases),  and  the  direction  of  the  long  ellipse^  axis 
does  not  as  a rule  coincide  with  the  direction  of  the  radiated  signal 
vector  E.  We  can  suggest  the  following  mechanisms  of  signal  depolar- 
ization during  its  vertical  propagation  in  the  glacier. 

The  first  mechanism  consists  in  anisotropy  of  the  electromagnetic 
properties  of  the  ice  crystals.  It  is  well  known  that  for  monocrystal 
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line  freshwater  ice  the  difference  in  the  refraction  coefficients  of 
the  extraordinary  i('Wo<j)|  and  ordinary  waves  at  the  10  KHz 

frequency  amounts  to  0.31  (fimc- n, ,0-5.00  — 4,69),  while  in  the  op- 

tical band  «neoc  — /i0o  = 1.314  - 1.309  = 0.005.  It  is  also  known  that 
the  crystals  in  the  glacier  may  have  an  ordered  arrangement,  spec- 
ifically: the  crystal  optical  axis  is  normal  to  the  shear  plane  of 

the  glacier  as  it  moves. 

Depolarization  of  the  propagating  wave  can  arise  in  those  cases 
when  the  wave  propagation  direction  does  not  coincide  with  the  crys- 
tal optical  axis  direction.  In  this  case  the  depolarization  is  ex- 
plained by  the  fact  that  the  radiated  wave  is  split  by  such  a medium 
into  quadrature  ordinary  and  extraordinary  waves,  propagating  with 
different  phase  velocities.  In  the  general  case  the  received  signal 
has  elliptic  polarization,  and  linear  polarization  is  possible  only 
if  the  difference  in  the  ordinary  and  extraordinary  wave  phases  at 
the  reception  point  is  a multiple  of  tt. 

Another  explanation  for  the  polarization  change  is  the  Faraday 
effect,  the  essence  of  which  is  that  in  certain  media  with  the  pre- 
sence of  a magnetic  field  coinciding  with  the  wave  propagation  direc- 
tion the  propagating  wave  remains  linearly  polarized  but  the  polar- 
ization phane  changes.  If  the  magnetic  field  direction  is  perpendic- 
ular to  the  wave  propagation  direction,  depolarization  also  shows  up. 
However  in  our  case  the  explanation  of  depolarization  by  the  Faraday 
effect  cannot  be  accepted  because  of  the  small  magnitude  of  this 
effect . 

There  is  still  another  mechanism  which  can  cause  depolarization  - 
the  stress  state  of  the  medium.  It  is  well  known  that  certain  sub- 
stances, specifically  glass,  in  the  presence  of  internal  stresses 
have  tensor  properties  which  are  not  associated  with  the  crystalline 
structure.  The  role  of  this  mechanism  in  the  glacier  is  quite  dif- 
ficult to  picture  and  evaluate  because  of  the  little  study  which  has 
been  made  ot  it. 
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Experiments  directed  toward  study  of  the  laws  governing  signal 
polarization  variation  in  vertical  radar  sounding  have  been  carried 
out  by  each  SAE,  beginning  with  the  15th.  Several  different  series 
of  measurement sjhave  been  made  and  have  shown  that  the  polarization 
diagrams  of  the  received  signals  are  different  for  different  radiating 
dipole  positions.  With  a certain  position  of  the  radiating  dipole 
(in  the  horizontal  plane)  we  observe  minimal  signal  depolarization  - 
the  reflected  signal  has  minimal  ellipticity  and  the  polarization 
plane  is  rotated  very  little  or  not  at  all  relative  to  the  radiated 
signal  polarization  plane. 

In  certain  experiments  the  problem  was  posed  of  establishing 
the  polarization  autocorrelation  radius  (at  least  with  respect  to 
received  signal  polarization  plane  rotation).  The  experiments 
showed  that  the  polarization  of  signals  with  different  time  lag  in 
the  general  case  is  different;  for  the  signal  with  constant  time 
lag  the  polarization  fluctuations  with  respect  to  polarization  plane 
rotation  angle  are  small.  For  example,  for  one  survey  line  at  a 
distance  of  20  meters  with  average  rotation  angle  35°  the  fluctuation 
amplitude  was  about  + 5 to  + 7° . We  note  that  the  degree  of  depolar- 
ization (increase  of  ellipticity  and  change  of  the  rotation  angle) 
is  larger,  the  greater  the  glacier  thickness. 

Alteration  of  signal  polarization  by  the  glacier  is  already  used 
in  practice.  By  using  crossed  dipoles  we  can  reduce  the  fluctuations 
of  the  signals  reflected  from  the  glacier  bed  and  attenuate  the  sig- 
nals reflected  from  nonhomogeneities  in  the  upper,  most  nonuniform 
layer  of  the  glacier.  The  use  of  crossed  dipoles  in  the  15th  and 
16th  SAE  gave  tangible  results.  Study  of  this  interesting  phenomenon 
is  continuing. 

Some  Results  of  Experimental  Studies  of  Overall  Attenuation. 

On  the  basis  of  the  above  discussion  we  can  conclude  that  the 
quantitative  estimate  (in  decibels)  of  overall  radio  signal  attenua- 
tion as  a function  of  the  radar  parameters  and  state  and  structure  of 
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Figure  8.  Theoretical  dependence  of  overall 
radio  signal  attenuation  on  glacier  thickness 
for  different  absorption  temperatures. 

1-  signal  attenuation;  2-  dB. 


the  glacier  can  be  made  using  the  formula: 


lg,CX+26l* 

v is  • ..*<J n «•*-.  ....  m.L  . ... IN ••  i 


where  the  quantity  40.7  dB,  defining  the  scattering,  focusing,  and 
depolarization  losses,  may  change  for  specific  conditions. 

The  results  of  Ne  calculation  using  (37)  are  shown  in  Figure  8 
Figure  9 shows  the  experimental  data  on  the  influence  of  surface 
thawing  on  signal  level  in  radar  sounding.  We  note  that  the  large 
signal  attenuation  increase  which  is  recorded  experimentally  upon 
melting  of  the  glacier  cannot  be  explained  by  the  reflection  coef- 
ficient increase  alone.  Thus,  even  a quarter-wave  water  layer  (e'M- 
on  a snow  surface  (e'H  «2j  has  reflection  coefficient  equal  to 

0.965,  which  leads  to  upper  boundary  reflection  losses  Np  of  approx 
imately  23  dB. 

The  overall  attenuation  is  in  general  a random  quantity  (see 
Section  6).  This  is  clearly  confirmed  by  the  reflected  signal 


Figure  9.  Influence  of  surface  thawing  on 
reflected  radio  signal  level,  a)  air  temperature 
in  the  observation  period;  b)  voltage  of  signal 
reflected  by  the  bed  (in  relative  units). 

1-  snow;  2-  rain;  3-  Jan.;  Feb. 


fluctuations.  Nevertheless  the  large  number  of  overall  attenuation 
measurements  makes  it  possible  to  perform  averaging  and  verify  the 
validity  of  the  adopted  theoretical  model  for  various  conditions 
(glacier  thickness,  glacier  average  annual  temperature,  equipment 
parameters).  The  contribution  of  the  individual  components  of  Ne 
can  be  determined  by  special  experiments. 


Unique  studies  were  made  in  the  14th  SAE  to  determine  the  over- 
all signal  attenuation  frequency  dependence.  The  experiments  were 
made  using  a specially  developed  radar  set  having  frequencies  of  60, 

70,  100,  140,  200,  and  400  MHz.  The  experiments  confirmed  the  valid- 
ity of  the  adopted  radio  signal  propagation  model  when  sounding 

glaciers.  In  accordance  with  this  model  we  can  calculate  the  geometric 

2 

j.  losses  using  the  1/h  law,  and  the  focusing  factor  and  the  losses  in 

reflection  from  the  glacier  bed  can  be  calculated  as  indicated  above, 
s 

The  losses  resulting  from  scattering  for  cold  glaciers  are  about  1 
f dB  and  Increase  to  20  - 30  dB  with  thawing  of  the  surface  layer. 
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To  calculate  the  absorption  magnitude  we  need  to  know  the  ab- 
sorption temperature,  which  can  very  from  0.75  to  0.9  of  the  annual 
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Figure  10.  Averaged  overall  attenuation  of  radio 
signals  of  different  frequency  along  a single 
measurement  line,  a)  60  MHz;  b)  70  MHz;  c)  100  MHz; 
d)  140  MHz;  e)  200  MHz;  f)  440  MHz;  g)  attenuation 
of  radio  signal  on  segment  A (to  enlarged  scale); 

1)  annual  average  air  temperature  along  line; 

2)  glacier  thickness:  horizontal  scale  is  distance 

south  of  Molodezh  Station. 
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average  air  temperature  near  the  surface,  which  does  not  lead  to 
large  errors  in  determining  the  overall  attenuation.  However  it  is 
better  to  solve  the  inverse  problem:  accounting  for  all  the  remain- 

ing components  of  Ne,  we  identify  from  the  measured  attenuation  the 
absorption  losses  and  determine  the  absorption  and  the  absorption 
temperature,  which  characterizes  the  temperature  regime  of  the 
glacier.  Repeated  experimental  studies  in  the  15th,  l6th,  and  17th 
SAE  confirmed  the  conclusions  drawn  above.  The  experimental  data 
were  used  to  plot  curves  of  overall  radio  signal  attenuation  for  six 
frequencies  (Figure  10),  which  show  weak  frequency  dependence  of  the 
overall  attenuation  for  the  Antarctic  glaciers. 


6.  STUDY  OF  REFLECTED  SIGNAL  FLUCTUATIONS . GLACIER  VELOCITY 
MEASUREMENT 


By  fluctuations  we  mean  a definite  part  of  those  changes  which 
the  signal  undergoes  as  it  passes  through  the  glacier  and  returns 
to  the  receiving  antenna.  Generally  speaking,  the  signal  amplitude, 
phase,  envelope  shape,  lag,  and  finally,  polarization  change. 


The  received  signal  envelope  shape  (including  "splitting"  of 
the  signal  into  two  or  more  signals)  and  also  the  signal  lag  relative 
to  the  time  of  emission  make  it  possible  to  define  stratified  gla- 
cier structure  or  bed  irregularities  with  the  degree  of  detail  pro- 
vided by  the  radar  range  and  angle  resolution  and  also  make  it 
possible  to  determine  the  depth  at  which  these  irregularities  lie. 
When  using  a radar  with  signal  detection  we  do  not  observe  phase 
changes . 

In  the  case  of  matched  polarization  of  the  transmitting  and  re- 
ceiving dipoles  the  signal  amplitude  varies  regularly  with  increase 
of  the  glacier  thickness  and  variation  of  its  temperature  and  random- 
ly with  respect  to  many  factors,  unambiguous  identif ication  of  which 
is  not  possible  at  the  present  time. 

In  this  section  we  shall  analyze  the  random  component  in  the 
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signal  amplitude  variation  which  is  determined  by  the  different 
influence  of  scattering  at  the  glacier  surface  and  at  its  lower 
boundary,  by  the  polarization  change,  by  the  difference  in  absorption 
in  micrononhomogeneities,  and  so  on.  Use  of  amplitude  samples  ob- 
tained in  accordance  with  the  Kotel'nikov  theorem  (r<7v:  amplitude 
recording  interval  T shorter  than  the  stationary  random  process  high- 
frequency  component  period)  detailed  statistical  analysis  of 

the  amplitude  fluctuations.  This  analysis  is  important  in  solving 
several  practical  problems.  For  example,  knowledge  of  the  autocor- 
relation radius  Ra  makes  it  possible  to  evaluate  the  effectiveness 
of  weak  signal  energy  accumulation  when  sounding  a glacier  from 
carrier  vehicles  traveling  comparatively  slowly  with  respect  to  the 
glacier  and  from  a fast-flying  airplane. 

Photographic  recording  of  the  reflected  radio  signal  fluctuations 
was  first  accomplished  by  the  l6th  SAE  (January,  1971).  Recording 
was  accomplished  at  three  frequencies  (60,  213,  440  MHz)  from  a 
slowly  moving  vehicle  on  segments  whose  lengths  in  the  various  exper- 
iments were  20  - 100  M.  The  recording  frequency  was  200  - 250  frames 
per  100  meters  of  survey  line.  Signal  amplitude  was  displayed  using 
type  A scan.  Since  the  electromagnetic  signal  alters  its  polarization 
as  it  propagates  in  the  glacier,  and  since  the  nature  of  the  sub- 
glacier bed  is  probably  different  in  the  direction  of  glacier  flow 
and  in  the  perpendicular  direction,  the  measurements  were  made  along 
lines  arranged  in  the  form  of  a cross.  Figure  11  shows  the  reflected 
signal  amplitude  as  a function  of  distance.  The  large  amount  of 
experimental  data  on  reflected  signal  fluctuation  was  machine  pro- 
cessed. For  this  purpose  the  filming  interval  was  first  normalized 
(the  vehicle  travels  at  different  speeds  in  different  series  of 
measurements).  The  computer  was  used  to  calculate  the  autocorrelation 
relations  for  the  signals  of  the  three  frequencies  for  all  the  survey 
lines  and  measurement  points  and  the  cross-correlation  relations  of 
the  fluctuating  signals  of  different  frequencies  along  the  same  survey 
line . 


Analysis  of  the  reduced  data  made  it  possible  to  determine  the 
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Figure  11.  Fluctuations  of  radio  signals  reflected 
by  bed.  a,  b)  at  6 0 MHz  for  two  mutually  perpen- 
dicular directions  with  Ra  = 5.2  and  5.1  m respect- 
ively; c)  at  213  MHz  with  Ra  = 1.6  m;  d)  at  440  MHz 
with  Ra  = 0.8  m.  1-  signal  voltage  amplitude;  2- 
distance  along  track. 


fluctuation  autocorrelation  radius.  For  the  60  MHz  frequency  the 
autocorrelation  radius  is  approximately  5.1  - 5.2  m,  for  the  213  and 
440  MHz  frequencies  it  is  1.6  and  0.2  m respectively.  It  was  also 
established  that  the  scattering  properties  of  the  glacier  bed  are  the 
same  for  different  measurement  survey  line  orientations;  the  fluc- 
tuation pattern  is  more  stable,  the  lower  the  frequency  and  the  great- 
er the  glacier  thickness. 

The  data  of  statistical  analysis  of  the  signal  fluctuations 
indicate  the  possibility  of  determining  glacier  surface  velocity 
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relative  to  the  fluctuation  pattern  fixed  with  the  bed.  For  this 
purpose  it  is  necessary  to  make  repeat  measurements  of  the  fluctua- 
tions on  the  survey  lines  with  some  time  interval.  If  we  then  super- 
pose the  fluctuation  diagrams  obtained  at  different  times  we  can 
determine  their  relative  shift,  which  is  proportional  to  the  glacier 
surface  dispacement  and  finally,  considering  the  time  interval  be- 
tween pictures,  we  can  determine. the  average  surface  velocity.  We 
assume  that  the  surface  velocity  measurement  accuracy  will  amount  to 
a fraction  of  the  forementioned  correlation  radii.  We  emphasize 
that  this  method  makes  it  possible  to  measure  the  glacier  velocity  in 
regions  far  from  exposed  rock  formations. 

7.  RADIOINTROSCOPY  OF  GLACIERS 

By  glacier  radiointroscopy  we  mean  study  of  the  internal  struc- 
ture of  glaciers  with  the  aid  of  radio  waves.  The  primary  elements  of 
glacier  microstructure  are  its  layering,  presence  of  foreign  inclu- 
sions in  its  interior,  continuity  of  the  cover,  and  structure  of 
the  lower  boundary  of  the  glacier,  which  Is  in  contact  with  the  bed. 

Study  of  Glacier  Stratified  Structure  and  Foreign  Inclusions. 

We  mentioned  previously  that  glaciers  have  three  primary  layers: 
snow,  firn,  and  ice.  The  radar  method  makes  it  possible  to  study  the 
nonhomogeneity  of  these  layers  and  detect  foreign  inclusions  in  the 
glacier,  for  example,  moraine-bearing  strata. 

Two  types  of  reflecting  nonhomogeneities  are  possible  in  the 
interior  of  the  glacier:  layers  with  different  glacier  material 

density  or  different  crystalline  structure,  and  the  ice  layers  which 
contain  rock  inclusions  - moraines.  Two  groups  of  reflected  signals 
correspond  (on  the  oscillograms  or  in  time)  to  these  nonhcmogeneit ies : 
one  is  near  the  sounding  pulse  and  the  other  is  near  the  pulse  re- 
flected by  the  bed.  We  note  that  reflected  signal  interpretation  is 
very  difficult;  a special  technique  and  specially  planned  experiments 
and  processing  of  the  results  are  required.  In  some  cases  the 
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reflected  signals  may  indicate  the  presence  of  moraine-bearing  layers 
near  the  bed.  For  example,  sounding  near  a borehole  in  the  Mirnyy 
Observatory  region  showed  that  the  glacier  thickness  here  is  120 
meters,  however  a reflected  signal  was  also  received  from  a depth 
of  about  80  m.  Later,  morainic  material  was  detected  during  drilling 
in  the  ice  at  a depth  of  75  m,  drilling  was  terminated  because  of  the 
presence  of  this  material. 


We  shall  present  some  other  examples. 


In  the  12th  SAE  during  simultaneous  sounding  of  the  glacier  at 
frequencies  of  100  and  ^ 4 0 MHz  layered  glacier  structure  was  re- 
corded at  a distance  of  30  - 80  meters  from  the  glacier  bed.  In 
1971  during  the  16th  SAE  a radar  operating  at  440  MHz  detected  on 
a survey  line  passing  to  the  south  of  Molodezh  Station  a nonhomo- 
geneity layer  about  5 km  long,  lying  at  a depth  of  about  1000  m, 
approximately  half  the  glacier  thickness. 


Radiosounding  of  the  glacier  from  an  airplane  on  a route  from 
Molodezh  Station  along  the  meridian  of  this  station  at  72°  S,  and 
then  after  a turn  to  the  Japanese  Syov  Station,  and  from  there  back 
again  to  the  Molodezh  Station  gave  unique  results.  Signals  whose 
time  lag  relative  to  sounding  pulse  is  greater  than  the  lag  of  the 
signals  reflected  from  another  interface  lying  below  the  glacier  bed 
were  recorded.  It  is  believed  that  this  interface  is  the  freezing 
boundary  of  the  bedrock  located  beneath  the  glacier.  It  should  be 
noted  here  that  correct  interpretation  of  the  signals  received  from 
nonhomogeneities  in  a glacier  is  aided  significantly  by  observations 
with  the  aid  of  boreholes,  observations  of  the  glacier  face  along 
the  coastline,  and  also  inverted  icebergs. 

Detection  of  Hidden  Crevasses. 

The  detection  of  crevasses  hidden  by  snow  bridges  is  a very 
important  problem,  since  its  successful  solution  makes  it  possible 
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to  ensure  safety  of  vehicle  movement  on  crevassed  glaciers.*  Evalua- 
tion of  ice  massif  continuity  on  areas  selected  for  takeoff  and 
landing  strip  construction  is  also  very  important.  Radar  crevasse 
finders  not  only  facilitate  considerably  the  selection  of  such  areas, 
but  also  make  possible  continuous  monitoring  of  their  condition.  An 
attempt  was  made  in  the  12th  SAE  to  detect  crevasses  from  an  airplane 
by  radar.  The  experiments  showed  that  the  signals  reflected  from 
crevasses  located  at  distances  up  to  four  kilometers  are  easily 
tracked  during  flight  at  300-meter  altitudes.  This  showed  tha  the 
signals  reflected  from  crevasses  can  be  detected  with  small  rf  wav 
glance  angles  (relative  to  the  glacier  surface). 

For  detecting  and  determining  the  location  and  orientation  of 
crevasses,  and  also  for  measuring  their  width,  it  is  best  to  use 
perspective  sounding  or  panoramic  radar  scan  from  a surface  vehicle 
or  an  airplane. 

The  overall  attenuation  of  the  signals  reflected  from  crevasses 
when  using  perspective  glacier  radiosounding  can  be  calculated  using 
the  general  scheme  presented  in  Section  5.  It  is  obvious  that  for 
detection  and  exact  determination  of  crevasse  location  and  dimen- 
sions, it  is  necessary  to  use  equipment  with  better  spatial  resolu- 
tion than  in  the  case  of  glacier  thickness  measurement,  for  example, 
with  radiation  pattern  width  10°  and  pulse  duration  10  ps.  However 
the  first  attempts  to  detect  crevasses  by  the  panoramic  radar  scan 
method  from  a surface  vehicle  were  made  with  the  aid  of  the  very 
simple  RV-10  equipment,  which  is  used  for  vertical  sounding. 

The  effective  reflecting  surface  area  of  glacier  crevasses  and 
the  backscatter  diagram  have  not  yet  been  studied.  Therefore  for 


* Crevasses  arise  basically  because  of  nonuniformity  of  glacier 
movement.  Their  depth  reaches  hundreds  of  meters  and  they  are 
covered  over  by  snow  bridges.  During  the  first  radiosounding  exper- 
iments the  tractor-sled  train  broke  through  into  a crevasse  22  meters 
deep . 


an  approximate  evaluation  we  can  consider  the  crevasse  to  be  a plane 
layer  with  rough  boundaries.  The  irradiated  crevasse  segment  area 
(in  m ) in  the  limits  of  a single  radar  spatial  resolution  element 
can  be  determined  from  geometric  considerations  using  the  following 
expression : 

Ci~*ma  l 1 1n  eirrma 

2 tstt  V H ~ - 

where  a = glance  angle  in  relation  to  the  crevasse  plane;  R = distance 
to  crevasse. 

Trepov  gave  the  parameters  of  the  equipment  used  for  the  experi- 
ments and  presented  the  necessary  calculations,  based  on  account  for 
the  radio  signal  attenuation  components  in  perspective  sounding  for 
the  case  with  antenna  located  8 meters  above  the  glacier  surface. 

When  using  radar  with  frequency  440  MHz  (G  = 2;  rBlu,  = 0.3  ys,  poten- 
tial 130  dB)  signals  reflected  from  crevasses  can  be  detected  at 
distances  over  100  meters  and  when  using  radar  with  frequency 
3000  MHz  (G  = 300,  t**b  = 10  ns)  and  the  same  potential  they  can  be 
detected  to  nearly  100  meters. 

The  problem  of  radar  detection  of  crevasses  is  complicated  by 
the  presence  of  noise  in  the  form  of  signals  reflected  from  irregular- 
ities on  the  glacier  surface  (sastrugi)  and  nonhomogeneities  in  the 
upper  part  of  the  glacier.  The  level  of  this  noise  is  determined  by 
the  glance  angles,  antenna  pattern  width,  degree  of  surface  irregular- 
ity, and  density  nonhomogeneity  of  the  glacier.  Experiments  are 
required  to  evaluate  the  noise  level. 

Experimental  studies  were  made  during  the  l4th  SAE  to  clarify 
the  possibilities  of  using  perspective  radar  sounding  for  crevasse 
detection  in  glaciers.  The  sounding  was  made  with  the  aid  of  the 
RV-10  radar.  The  indicator  was  an  SI-20  oscillograph.  The  receiving 
and  transmitting  antennas  were  mounted  on  a mast  on  the  roof  of  the 
mobile  laboratory  at  a height  of  7 - 8 meters  above  the  glacier 
surface.  The  studies  were  made  near  Molodezh  Station  on  the  airfield 
and  the  glaciological  test  range,  and  also  in  the  crevasse  zone  near 
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Slbiryachka  Bay  (March,  1969).  The  oscillograms  show  a significant 
reflected  signal  level  at  distances  0 - 300  M when  sounding  in  tne 
crevasse  zone  and  practically  no  signal  in  regions  with  monolithic 
ice  (airfield).  It  was  noted  that  the  reflected  signal  level  is 
considerably  higher  with  the  antenna  oriented  perpendicular  to  the 
crevasses  than  with  orientation  along  the  crevasses  (Figure  12). 

Thus  the  radar  signals  reflected  from  crevasses  are  observed 
reliably  on  the  perspective  sounding  radar  screen.  However,  because 
of  the  low  spatial  resolution  of  the  equipment  used,  it  is  rarely 
possible  to  establish  unique  correspondence  between  the  radar  signals 
and  crevasse  location  in  the  glacier.  The  use  of  equipment  with 
higher  spatial  resolution  will  obviously  make  it  possible  to  solve 
this  problem.  Modified  equipment  was  tried  out  in  the  l8th  SAE 
(1973).  Experiments  on  perspective  crevasse  sounding  were  made  at 
440  MHx  frequency  with  use  of  an  antenna  array  consisting  of  eight 
radiators  with  a flat  reflector.  This  equipment  was  developed  by 
the  radiophysics  division  of  the  Arctic  and  Antarctic  Scientific 
Research  Institute  (AANII)  in  1972. 

Study  of  Glacier  Bed  Relief. 

The  relief  of  the  lower  boundary  of  the  glacier  is  completely 
identical  to  the  relief  of  the  underlying  bedrock.  Therefore  study 
of  the  structure  of  the  relief  of  the  lower,  bottom  surface  of  the 
glacier  reduces  to  study  of  its  sub-ice  relief. 

Study  of  the  sub-ice  relief  was  initiated  by  the  10th  SAE  and 
is  continuing  at  the  present  time.  The  results  of  the  studies  are 
used  by  the  Division  of  Aerial  Photography  and  Cartography  of 
Soyuzmorniiproekta  and  by  the  Division  of  Geography  of  Polar  Countries 
of  AANII  for  mapping  the  sub-ice  relief.  Recording  of  the  sub-ice 
relief  includes  measurement  of  airplane  altitude  relative  to  sea 
level  , radar  measurements  of  flight  altitude  above  the  glacier 

i (//p-t  » and  radar  measurements  of  glacier  thickness  (h).  These  data 


are  used  to  calculate  the  position  of  the  glacier  surface  (A„  = Hn  — H^) 
and  the  glacier  bed  ( A. 7 *//«--  h)  relative  to  sea  level. 


Figure  12.  Oscillograms  of  perspective  glacier 
radiosounding,  a)  sounding  in  crevasse  region 
near  Sibiryachk  Bay,  antenna  oriented  along 
crevasses;  b)  same,  antenna  orient  d across 
crevasses;  c)  sounding  in  region  without  crevasses 
(airfield,  Molodezh  Station).  The  signal  reflected 
from  the  glacier  bed  is  seen  on  the  right  side  of 
the  oscillograms. 
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The  radiosounding  survey  line  is  "referenced"  to  the  map  on  the  basis 
of  airplane  navigational  instrument  indications  and  data  of  aerial 
photography  of  the  point  of  airplane  approach  to  the  sounding  line 
from  the  coastline  and  the  point  of  departure  from  the  pass. 

The  degree  of  detail  of  the  sub-ice  relief  record  is  very  high 
(frame-by-frame  recording  of  the  radar  screen  is  accomplished  approx- 
imately every  kilometer).  This  detail  makes  possible  more  precise 
analysis  of  certain  characteristic  features  of  the  sub-ice  relief  and 
permits  more  valid  generalization  of  the  data  than  is  possible  when 
using  a large-scale  map  prepared  on  the  basis  of  radiobarometric 
surveying.  At  the  present  time  radar  surveying  has  been  performed 
on  routes  whose  overall  length  is  more  than  100,000  KM.  The  routes 
cover  an  area  of  about  one  million  square  kilometers. 

The  most  characteristic  and  prominent  feature  of  the  sub-ice 
relief  of  the  part  of  Antarctica  which  we  have  studied  is  its  irreg- 
ularity; we  frequently  encounter  depressions  and  peaks  whose  base 
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width  is  equal  to  or  somewhat  greater  than  their  vertical  dimension. 

This  contradicts  the  concepts  on  the  basis  of  which  the  existing  map 
of  the  sub-ice  relief  of  Antarctica  was  drawn  [2]. 

. 

The  radiobarometric  survey  data  led  to  the  detection  of  a "sub- 
ice" strait,  extending  from  Edward  VIII  Bay  to  Amudsen  Bay,  a large 
mountain  massif  to  the  south  of  this  "strait"  between  68  and  69.5°  S, 
a sub-ice  mountain  lowland  adjacent  to  the  Syov  Station,  a second 
large  mountain  massif  to  the  south  of  Syov  Station,  extending  south- 
ward from  72°  S latitude. 

In  the  17th  SAE  (1972)  a detailed  area-type  survey  of  the  sub- 
ice  relief  was  conducted  in  a special  tractor-sled  expedition.  Two 
test  ranges  were  selected  for  this  pumose  on  the  Mirnyi-km  153  route 
in  the  direction  of  Pioner  Station.  These  test  ranges  had  areas  of 

p 

5 x 5 KM  and  were  located  at  km  57  and  km  153  of  the  indicated  route. 
Radiosounding  was  accomplished  by  surveying  the  test  ranges  along 
mutually  perpendicular  lines  at  intervals  of  500  meters. 

The  survey  made  on  these  test  ranges  not  only  made  it  possible  to 
obtain  detailed  two-dimensional  maps  of  the  sub-ice  relief,  but  also 
made  it  possible  to  "tie  in"  to  the  bed  relief  the  test  range  points 
marked  by  stakes. 

Figures  13  and  14  illustrate  the  work  done  in  studying  the 
spatial  structure  of  the  lower  boundary  of  the  glacier.  Repeat 
surveys  on  these  test  ranges  will  make  it  possible  to  evaluate  glacier 
displacement  relative  to  the  "bottom"  reference  points. 

In  conclusion  we  note  that  radiosounding  makes  it  possible  to 
obtain  the  very  detailed  information  which  is  necessary  for  map 
•>  preparation  and  geologic  study  of  the  bedrock. 
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8.  RADIOSOUNDING  OF  TEMPERATURE  AND  MOUNTAIN  GLACIERS 


Characteristics  of  Glacier  Structure  and  Properties. 

We  have  discussed  the  theoretical  and  experimental  studies 
directed  toward  the  development  and  introduction  of  radar  methods 
for  extensive  study  of  the  cold  glaciers  of  Antarctica.  Radio- 
sounding of  the  Arctic  (temperate)  and  mountain  glaciers  is  also 
of  great  practical  interest.  The  radar  methods  are  without  question 
the  most  precise  and  fastest  and  can  be  used  to  evaluate  the  water 
resources  of  these  glaciers  and  study  their  physical  characteristics. 

The  temperate  and  mountain  glaciers  differ  markedly  from  the 
Antarctic  glaciers  in  their  properties:  their  average  interior 

temperature  is  higher  and  sometimes  close  to  0°C,  on  the  glacier 
surface  there  is  nearly  always  a layer  of  wet  snow  or  water,  within 
the  glacier  there  are  layers  of  different  density  corresponding 
to  the  annual  cycle  of  glacier  mass  accumulation  and  layers  of  foreign 
matter  due  to  rock  formation  weathering.  The  specific  properties 
of  the  mountain  glaciers  make  the  radar  method  in  the  form  in  which 
it  is  used  in  Antarctica  ineffective  for  measurement  of  the  thickness 
of  these  glaciers. 

It  is  well  known  that  electromagnetic  wave  dielectric  losses 
increase  markedly  with  increase  of  the  ice  temperatures  to  values 
close  to  0°C.  Table  7 presents  the  energy  absorption  of  a plane 
wave  as  it  propagates  over  a distance  of  two  kilometers  (glacier 
thickness  one  kilometer).  The  layer  of  wet  snow  or  water  on  the 
surface  of  the  thawing  glacier  leads  to  significant  reflection  from 
the  air-ice  interface  in  comparison  with  the  dry  glacier  case. 

While  at  the  air-ice  interface  the  reflection  coefficient  is  0.28, 
with  a quarter-wave  water  layer  present  on  the  surface  of  the  ice 
the  coefficient  is  0.92  (without  account  for  the  dielectric  losses). 
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TABLE  7 


Frequency,  MHz  .... 

...  440 

. 200 

100 

60 

30 

10 

Absorption,  dB 

at  -T’  C 

, . . 102.4 

97.8 

91.2 

87 

81 

70 

- 

at  - 20°  C About  30 


Increase  of  the  attenuation  of  the  440  MHz  frequency  radar  sig- 
nal upon  glacier  thawing  was  demonstrated  experimentally  in  the  11th 
and  12th  SAE:  with  a wet  snow  layer  10  - 15  CM  thick  present  the 

attenuation  was  20  - 30  dB.  The  signal  attenuation  because  of  re- 
flection from  layers  differing  in  electrical  properties  (different 
ice  density,  rock  inclusions)  depends  on  the  glacier  structure. 
According  to  measurements  made  in  the  12th  SAE  (1967)  at  this  same 
frequency  the  attenuation  is  only  a few  decibels  for  the  Antarctic 
glaciers.  For  the  more  nonhomogeneous  mountain  glaciers  the  attenua- 
tion because  of  internal  reflections  is  considerably  higher.  With 
the  presence  of  periodic  layers  commensurate  in  thickness  with  the 
wavelength  the  attenuation  may  reach  tens  of  decibels. 


From  comparison  of  the  overall  radio  signal  attenuation  in 
mountain  and  Antarctic  glaciers,  we  can  conclude  that  the  radar 
with  carrier  frequency  440  MHz  and  receiving  channel  sensitiviey 
130  dB  relative  to  the  peak  transmitter  radiated  power  can  measure 
mountain  glacier  thickness  reliably  to  80  - 100  M.  Considering 
that  mountain  glacier  thickness  reaches  300  meters  or  more,  lower 
frequencies  (10  - 100  MHz)  should  be  used. 


Another  characteristic  of  mountain  glaciers  is  the  marked  ir- 
t regularity  of  their  surface.  In  this  connection  reflections  from 

their  surface  irregularities  are  possible  when  radiosounding  such 
glaciers.  These  reflections  made  Identification  of  the  useful  signal 

k 

difficult.  Highly  directive  antennas  must  be  used  to  reduce  the 
interfering  reflections. 
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Calculation  shows  that  the  radar  set  with  receiver  sensitivity 
150  dB  relative  to  the  peak  transmitter  power,  carrier  frequency 
60  MHz,  and  10  dB  gain  antenna  can  provide  mountain  glacier  thickness 
measurement  to  within  300  - 400  M. 

Radiosounding  of  Soviet  Arctic  Glaciers. 

In  1968  the  Sever-20  High-Latitude  Expedition  conducted  the  first 
experiments  in  the  U.S.S.R.  on  radiosounding  of  the  Severnaya  Zemlya 
Archipelago  glaciers  [11].  English  and  American  scientists  made 
similar  studies  in  1966,  performing  airborne  radar  measurements  of 
glacier  thickness  on  Ellsmere  Island. 

The  arctic  glaciers  have  comparatively  small  dimensions  and  many 
landmarks  are  available  in  the  periglacial  zone.  Therefore  aerial 
photography  makes  it  possible  to  obtain  precise  tie-in  of  the  sounding 
survey  lines.  The  Severnaya  Zemlya  Archipelago  glacier  thickness  does 
not  exceed  1000  meters.  The  regional  glaciological  conditions  and 
the  meterological  characteristics  of  certain  regions  of  the  Arctic 
require  selection  of  the  most  suitable  time  of  the  year  for  sounding, 
when  meltwater  is  not  present  in  the  upper  glacier  layers,  the 
weather  is  sunny,  the  downslope  winds  are  minimal,  and  consequently 
low-altitude  flights  above  the  glacier  are  safest.  The  limited  time 
with  favorable  weather  conditions  in  the  period  in  which  the  fore- 
mentioned  experiments  were  conducted  shortened  the  duration  of  the 
studies  and  therefore  not  all  the  most  characteristic  Severnaya 
Zemlya  Archipelago  glaciers  were  studied  equally  thoroughly.  Figure 
15  shows  a schematic  of  flights  made  for  glacier  radiosounding.  We 
see  from  the  figure  that  radar  thickness  measurements  were  made  on 
the  glaciers  of  Schmidt,  Komsomolets,  Pioner,  Oktyabr'skoe  Revolyutsii 
Islands,  forming  the  Severnaya  Zemlya  Archipelago,  and  on  the  glacier 
of  Ushakov  Island,  located  in  the  Kara  Sea  to  the  northwest  of  this 
archipelago. 

A radar  with  the  following  characteristics  was  used  for  the 
radiosounding : 
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Figure  15.  Diagram  of  flights  made  in  radiosounding 
Severnaya  Zemlya  Archipelago  glaciers.  1-  Schmidt  G.; 
2-  Academy  of  Sciences  G. ; 3-  Rusakov  G. ; 4-  Pioner  G.  ; 
5-  Al'banov  G.;  6-  Ushanov  G. 


Carrier  frequency f = 100  MHz 

Pulse  power F„„n  = 5-7  W 

Pulse  duration  • • • • • fitMn  = 0*3  U S 

Receiver  sensitivity Pm  = 5 • 10~^  W 

Separate  "double  square"  antennas,  located  on  the  fuselage  of  an 
Li-2  airplane  in  the  wing  region  and  spaced  two  meters  from  one 
another,  were  used  to  radiate  and  receive  the  signals.  The  flight 
altitude  above  the  glacier  surface  was  measured  by  this  same  radar. 
Barometric  measurements  were  made  in  order  to  reference  the  air- 
glacier,  air-rock  formation,  and  glacier  glacier-bed  interfaces  to 
sea  level.  The  baromatric  altimeter  indications  were  "referenced" 


to  sea  level  in  the  regions  of  approach  to  the  coastline  and  departure 
from  the  coastline  based  on  the  airplane  radio-altimeter  indications. 
The  measurements  were  recorded  by  an  automatic  photocamera. 

It  seems  to  us  that  the  best  time  of  the  year  for  radar  sounding 
of  Arctic 'glaciers  is  spring  (March,  April,  early  May).  In  this 
period,  glacier  thawing  is  still  not  marked  and  electromagnetic  wave 
energy  absorption  is  not  significant.  Another  important  factor  are 
the  long  periods  of  sunny  weather  — they  ensure  flight  safety  at  low 
altitude  above  the  glacier.  The  studies  undertaken  by  the  Sever-20 
expedition  made  it  possible  to  obtain  important  information  on  thick- 
nesses of  the  glaciers  studied  with  error  no  more  than  -7,  + 3%, 
define  more  exactly  the  glacier  bedrock  location  in  relation  to  sea 
level,  and  evaluate  the  relief  of  this  bedrock. 

Radar  signal  attenuation  calculations  showed  that  during  conduct 
of  the  measurements  the  electromagnetic  signal  losses  in  the  glacier 
due  to  absorption  do  not  exceed  those  for  the  continental  Antarctic 
glaciers  in  the  10  - 15  kilometer  zone  along  the  coast.  All  the 
glacier  thickness  data  are  summarized  in  the  tables  of  [11]. 

9.  RADIOSOUNDING  OP  SEA  AND  FRESHWATER  ICE 

By  radiosounding  of  sea  and  freshwater  ice,  we  mean  remote 
determination  of  the  thickness  and  study  of  the  structure  of  the  ice 
covers  of  various  water  bodies  with  the  aid  of  active  radar  methods. 

We  believe  that  this  problem  was  first  formulated  in  [7].  However 
the  practical  use  of  the  results  of  that  study,  obtained  for  sim- 
plified models  of  freshwater  ice,  was  limited. 

The  importance  of  developing  a reliable  remote  method  and  radar 
equipment  for  sounding  floating  ice  covers  is  obvious.  Expansion  of 
the  commercial  transport  volume  in  the  Arctic,  study  and  practical 
use  of  the  Arctic  sea  shelf,  navigation  in  the  northern  fiver  estuaries 
scientific  prognosis  of  ice  conditions  — all  these  are  very  important 
economic  problems.  Their  successful  solution  requires  methods  for 
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fast  collection  of  data  on  the  ice  cover  thickness  over  large  areas. 
The  availability  of  such  hardware  aboard  the  airplanes  used  for  ice 
scouting  and  navigation  support,  and  also  aboard  the  helicopters  of 
the  icebreakers,  would  broaden  markedly  the  information  on  ice  cover 
status  in  their  operating  zone  and  would  facilitate  the  selection 
of  ship  navigation  and  maneuvering  tactics. 


It  is  well  known  that  the  development  of  radar  for  sounding 
the  ice  covers  of  different  water  bodies  is  a very  complex  physical 
and  engineering  problem.  Its  complexity  is  determined  primarily 
by  the  absence  of  generalizing  electromagnetic  characteristics  for 
sea  and  freshwater  ices  in  the  meter  to  decentimeter  wavelength  band. 
The  ice  covers  of  water  bodies,  particularly  the  Arctic  seas,  have 
very  complex  structure:  they  include  pack,  perennial,  and  young  ice 

of  differing  thickness.  All  these  different  ice  varities  have  dif- 
ferent electromagnetic  characteristics,  varying  markedly  across  their 
section  as  a function  of  ice  age,  thickness,  and  temperature. 
Therefore  statistics  on  ice  cover  geometric  (dimensional)  parameters 
and  a large  amount  of  data  characterizing  the  variability  of  their 
temperature  and  salinity  conditions  are  necessary  for  optimization  of 
the  technical  solutions. 


It  is  obvious  that  particularly  important  in  this  problem  are 
reliable  data  on  floating  ice  cover  electromagnetic  characteristics 
(e'  and  tg6)  and  the  physics  of  the  so-called  matching  layer,  at  the 
expense  of  which  ice  incrementation  takes  place.  This  layer  is 
extremely  nonhomogeneous , both  as  a dielectric  and  as  material  of  the 
phase  change  zone. 
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There  are  also  definite  complexities  in  the  technical  aspects. 
The  radar  must  have  high  time  resolution,  high  potential,  and  narrow 
antenna  pattern.  In  the  last  three  to  five  years  extensive  studies 
have  been  conducted  in  several  countries  and  particularly  in  the 
U.S.S.R.  of  floating  ice  cover  electromagnetic  characteristics, 
which  have  advanced  considerably  the  solution  of  the  problem  of 
developing  the  radar  method  for  sounding  floating  and  particularly 
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Arctic  drifting  ice. 

Electromagnetic  Characteristics  of  Sea  and  Freshwater  Ice. 

A quite  wide  frequency  band  is  applicable  for  sounding  sea  ice 
by  active"  radar  methods;  from  30  to  1000  MHz.  In  [9]  the  radar  range 
equation  was  used  to  calculate  the  overall  attenuation  of  electro- 
magnetic waves  during  propagation  in  sea  ice.  The  calculation  results 
are  shown  in  Figure  16.  1$  the  calculations  it  was  assumed  that  the 

Fresnel  coefficient  of  reflection  from  the  lower  boundary  is  equal  to 
one,  the  ice  layer  is  plane-parallel,  and  wave  incidence  is  normal. 
These  data  show  that  a radar  potential  of  140  - 160  dB  is  required 
for  sounding  sea  ice.  Full-scale  experimental  data  on  the  magnitudes 
of  the  refraction  index  n and  electromagnetic  energy  absorptivity  /V„ 
(dB/m),  which  determine  the  wave  phase  velocity  and  the  required  radar 
energy  potential,  are  necessary  for  optimal  solution  of  the  problem. 


At  the  present  time  results  are  already  available  of  laboratory 
and  full-scale  measurements  of  these  characteristics.  We  consider 
the  most  reliable  measurements  to  be  those  made  by  Khokhlov  on  the 
SP-13  drifting  station  [12]  and  those  made  by  Tripol'nikov  in  1972-73 
on  the  SP-21  drifting  station.  Khokholov  measured  e'  and  tg<$  of 
drifting  ice  in  the  100  Hz  to  1 MHz  band  and  showed  their  functional 
dependence  in  the  anomalous  dispersion  region  on  ice  salinity, 
structure,  and  temperature.  Tripol'nikov  studied  the  "quiet"* 
frequency  band:  30  - 400  MHz. 

Laboratory  and  full-scale  studies  of  e'  and  tg6  in  the  10  Hz 
to  100  MHz  band  are  discussed  in  [28-30].  These  data  do  not  cover 
the  entire  frequency  band  which  we  believe  to  be  promising  for  radar 
study  of  sea  ice.  It  should  be  emphasized  that  sea  ice  is  a strongly 
absorbing,  stratified  nonhomo geneous  dielectric  which  changes  its 


* By  "quiet"  we  mean  the  frequency  band  where 
approximately  constant  value. 


retains  an 
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properties  sharply  under  the 
influence  of  various  inputs  and 
is  very  difficult  to  model 
under  laboratory  conditions. 

In  addition  the  technique  of 
these  studies  is  based  on 
measuring  the  parameters  of 
concentrated  chains.  For  the 
meter-centimeter  wavelength 
band  the  estimates  of  e'  and 
tg<$  are  best  made  by  modeling 
the  wave  propagation  process 
in  sea  ice.  In  [13]  an  attempt 
was  made  to  measure  the  absorp- 
tion under  full-scale  conditions 
at  frequencies  of  100,  *440, 
and  1000  MHz.  These  results 
do  not  agree  with  the  data  of 
[28].  In  the  following  we 
shall  discuss  the  results  of 
full-scale  measurements  made 
by  Tripol 1 nikov. 


Figure  16.  Calculated  overall 
radio  signal  attenuation  in  sea 
ice  for  various  absorptivities . 

1-  antennas  at  one-meter  height; 

2-  antennas  at  100-meter  height; 

3-  overall  attenuation;  4-  ice 
thickness;  5-  dB;  6-  dB/m. 


The  data  of  electromagnetic  wave  propagation  velocity  and  ab- 
sorptivity measurements  in  the  30  - 400  MHz  band  in  sea  ice  were  ob- 
tained by  pulse  "sounding"  of  coaxial  long  lines  having  sea  ice  as 
the  dielectric  under  study.  In  such  a line  in  the  absence  of  condi- 
tions for  resonant  energy  exchange  between  the  field  components  E 
and  H the  signal  attenuation  is  due  basically  to  losses  in  the  di- 
electric. In  this  case  the  measured  specific  attenuation  can  be  taken 
as  the  maximal  estimate  of  the  specific  absorption.  It  will  be  clear 
from  the  quantitative  results  presented  below  that  when  using  sea  ice 
as  the  dielectric  in  the  line  the  specific  attenuation  practically 
coincides  with  the  specific  absorption,  and  the  wave  velocity  in  the 
line  coincides  with  the  velooity  in  the  medium.  Coaxial  lines  of 
length  from  two  to  eight  meters  were  used  for  the  measurements. 
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The  measurement  technique  was  developed  and  the  line  parameters 
were  optimized  under  laboratory  conditions.  Observations  were  made 
of  the  passage  of  radio'  signals  of  duration  from  0.05  to  0.5  usee 
with  frequencies  30,  60,  100,  200,  400  MHz  through  lines  filled  with 
ice  of  different  salinity,  and  also  video  pulses  of  duration  0.02 
usee  with- repetition  frequency  10  kHz.  The  wave  propagation  velocities 
were  measured  on  the  basis  of  the  delay  time  of  the  signal  passing 
through  the  line.  The  specific  attenuation  measurements  were  made 
by  the  variable  length  line  method  and  on  constant-length  lines  with 
account  for  the  losses  in  reflection  from  the  ends  of  the  line.  The 
measurements  were  made  at  temperatures  from  -40  to  -5°C.  Pulse  front 
distortion  was  observed  with  salinity  17%  and  temperature  above  the 
eutectic  for  NaCl.  But  since  the  group  velocities  coincide  for  the 
rf  and  video  pulses  this  effect  can  be  explained  by  the  frequency 
dependence  of  tg6 . 

It  was  also  noted  that  with  temperature  increase  above  -22°C, 
there  is  reduction  of  the  velocity,  and  the  more  marked  the  higher 
the  sample  salinity.  This  phenomenon  can  be  explained  by  the  appear- 
ance of  a liquid  phase  in  the  medium  under  study  — brine  droplets, 
whose  concentration  in  the  total  ice  volume  increases  with  increase 
of  the  temperature  and  salinity.  In  this  case  the  sounded  medium  is 
a mixture  of  ice  and  water  with  some  effective  e*.  This  effect  is 
of  fundamental  importance  for  correct  interpretation  of  sea  ice  radar 
sounding.  For  example,  the  annual  average  temperature  of  drifting 
sea  ice  at  a depth  equal  to  half  its  thickness  lies  approximately  in 
the  range  from  -5  to  -8°C,  the  salinity  of  year-old  ice  at  this  same 
level  is  about  4-5?.  For  such  ice  the  sounding  pulse  velocity  in 
the  ice  can  be  taken  as  140  m/usec.  In  perennial  pack  ice  the  salin- 
ity at  the  0.5  thickness  level  is  about  0.5%;  in  such  ice  the  elec- 
tromagnetic wave  velocity  is  170  m/usec.  The  series  of  laboratory 
measurements  made  it  possible  to  evaluate  the  specific  absorption  of 
electromagnetic  waves  in  ice  of  different  salinity  and  temperature. 

It  was  found  that  ice  phase  composition  variation  as  a function  of 
thickness  can  have  significant  influence  on  sounding  result  interpre- 
tation accuracy.  Therefore  it  becomes  necessary  to  obtain  the 


averaged  electromagnetic  characteristics  of  sea  ice  of  different  age 
and  ice-formation  conditions. 

We  have  mentioned  previously  the  need  for  studying  the  electri- 
cal properties  and  in  general  the  physics  of  the  matching  layer. 

In  the  present  case  it  is  very  important  to  explain  the  presence 
of  electromagnetic  contrast  at  the  crystallization  boundary  of  sea 
water.  This  question  is  controversial,  and  for  its  concretization 
we  shall  present  some  information  on  sea  water  crystallization  process 

The  first  ice  crystals  form  upon  cooling  sea  water  to  -1.9 °C. 
Further  temperature  reduction  facilitates  the  crystallization  process, 
as  a consequence  of  which  the  salt  concentration  in  the  solution 
increases.  The  concentrated  solution  (liquid  phase)  which  remains 
in  the  crystallization  process  fills  either  local  volumes  or  inter- 
crystalline layers,  forming  communicating  channels.  At  low  tempera- 
tures the  local  volumes  dominate  and  at  higher  temperatures  — the 
channels.  At  certain  temperatures  the  solution  can  become  saturated 
for  some  component  and  then  this  component  precipitates.  Thus  KC1 
precipirates  at  -11.1°C,  NaCl  at  -22.6°C.  CaCl ^ has  the  lowest 
eutectic  temperature,  it  precipitates  at  -55°C.  It  follows  from  this 
that  natural  sea  ice  cannot  contain  a liquid  phase  in  its  structure. 
The  crystallization  process  leads  to  capture  by  the  ice  crystal 
lattice  of  the  solution  ions.  And  although  their  concentration  is 
extremely  low  the  influence  of  impurities  on  the  electrical  charac- 
teristics of  crystallized  ice  is  very  significant.  The  influence  of 
the  liquid  phase  on  the  physical  properties  of  sea  ice  is  significant; 
at  not  very  low  temperatures  brine  occupies  a considerable  fraction 
of  the  total  volume.  The  temperatures  of  the  lower  layers  of  the 
formed  ice  cover  of  thickness  about  0.5  meters,  covered  with  a snow 
layer,  are  close  to  the  crystallization  temperature;  they  are  very 
strongly,  nearly  half,  saturated  with  concentrated  brine.  Therefore 
it  is  very  important  not  only  to  explain  the  presence  of  electro- 
magnetic contrast  of  this  unique  boundary,  but  also  to  estimate  its 
quantitative  characteristics  for  the  meter  and  decimeter  wavelength 
band.  We  have  not  found  such  information  in  the  literature. 
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Tripol'nikov  carried  out  a major  series  of  studies  on  the  SP-21 
drifting  station,  using  the  long-line  measurement  method,  This 
technique  made  it  possible  to  determine  electromagnetic  signal  atten- 
uation, wave  phase  velocity  in  the  ice,  and  the  presence  of  electro- 
magnetic contrast  between  the  ice  and  water.  The  experimental  study 
of  these  important  parameters  was  made  as  follows. 

Two  10  mm-diam  clear-through  parallel  holes  were  melted  in  sea 
ice  10  cm  from  one  another  by  a thermal  needle.  Metal  rods  5 mm  in 
diameter  with  length  exceeding  the  ice  thickness  were  inserted  into 
the  holes.  Their  lower  ends  extended  beyond  the  ice-water  inter- 
face. The  metal  rods  were  frozen  into  the  ice,  forming  a two-conduc- 
tor "ice"  line  with  wave  impedance  about  200  Ohms,  shorted  at  the 
lower  boundary  of  the  ice  being  studied  by  the  sea  water.  (It  was 
found  that  the  technology  used  to  make  the  holes  does  not  alter 
appreciably  the  ice  structure  between  "ice"  line  rods.)  Then  sounding 
of  this  line  was  carried  out  by  video  and  rf  pulses  with  the  aid  of 
a system  having  a delay  line  matched  with  the  "ice"  line.  The  pulses 
were  recorded  by  combined  exposure:  the  pulse  reflected  from  the 

short  circuit  in  the  line  made  by  the  operator  at  the  ice-air  inter- 
face was  exposed  first,  then  the  short  circuit  in  the  air  was  removed 
and  the  pulses  were  exposed  with  reflection  from  the  short  circuit 
at  the  lower  end  of  the  "ice"  line.  Additional  gain  was  not  inserted 
into  the  circuit  in  this  process. 

Comparing  the  signals  reflected  from  the  artificial  short 
circuit  at  the  air-ice  interface  and  from  the  short  circuit  formed 
by  the  sea  water,  we  can  evaluate  the  overall  losses  in  the  "ice" 
line.  These  losses  are  due  to  absorption  in  the  dielectric  (natural 
sea  ice)  and  by  reflection  at  the  air-ice  and  ice-water  interfaces. 

The  latter  losses  can  be  considered  a measure  of  the  electromagnetic 
contrast.  The  measurements  were  made  on  perennial  pack  ice  of 
thickness  up  to  10  meters,  on  two-year-old  pack  ice  of  thickness  from 
two  to  three  meters,  and  on  young  ice  one  meter  thick.  These  ice 
types  are  characterized  by  different  temperature  and  salinity  dis- 
tribution through  the  thickness.  The  experimental  values  of  att  - 
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tion  as  a function  of  ice  thickness  are  shown  in  Figure  17.  Here 
the  dashed  curves  are  the  theoretical  curves  corresponding  to  part- 
icular specific  absorption  values. 

In  spite  of  the  considerable  variety  of  the  general  physical 
characteristics  of  the  ices  studied,  the  results  of  sounding  by 
25-ns-duration  video  pulses  show  that  in  most  cases  the  specific 
absorption  lies  in  the  range  from  0.9  to  2.0  d3/m.  The  reflectivity 
can  be  taken  equal  to  one,  otherwise  the  attenuation  curves  would 
cross  the  abscissa  axis  above  1.0.  This  is  also  valid  for  200-MHz- 
frequency  rf  pulses  and  the  absorptivity  is  4-5  dB/m.  Confirmation 
of  the  existence  of  marked  electromagnetic  contrast  at  the  ice-water 
interface  was  also  obtained  in  sounding  lines  by  rf  pulses  at  4^"o 
MHz.  However  quantitative  evaluation  at  440  MHz  with  the  line  variant 
is  difficult  — reduction  of  the  distance  between  the  "ice"  line 
conductors  requires  a more  compelx  thermal  drilling  technique.  The 
results  shown  here  do  not  indicate  precisely  what  is  the  reflecting 
boundary  — the  sea  water  itself  or  the  upper  boundary  of  the  layer  in 
which  high-concentration  brine  is  present.  Radar  sounding  of  an 
"ice  island"  was  performed  to  obtain  an  answer  to  this  question 
(SP-19  drifting  station).  Its  thickness  was  28  - 3^  M.  Salinity 
samples  were  taken  from  various  depths  of  this  "island"  and  tempera- 
ture measurements  were  made.  The  salinity  and  temperature  values 
are-  shown  in  Table  8. 

Since  the  temperature  at  depths  of  17  meters  and  below  exceeds 
the  eutectic  temperature  for  all  the  seawater  salts,  at  these  depths 
the  liquid  phase  is  present  in  the  body  of  the  flow.  Radiosounding 
of  this  flow  was  carried  out  at  200  MHz  by  50  ns  duration  pulses. 

The  experiments  showed  that  the  salinity  boundary  located  at  the 
17  meter  depth  yields  a strong  reflected  pulse.  The  salinity  jump 
at  the  19  - 20  meter  level  also  reflects  a considerable  portion  of 
the  sounding  pulse  energy. 

Regular  reflection  was  noted  from  the  ice-water  interface.  This 
series  of  field  measurements  showed  that  electromagnetic  contrast 
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Figure  17.  Experimental  radio  signal  specific 
attenuation  as  function  of  ice  thickness. 

1-  25  ns  duration  video  pulses;  2-  50  ns  duration 
rf  pulses;  3-  specific  attenuation;  4-  ice  thickness; 
5-  dB/m. 


arises  in  ice  layers  characterized  by  variation  of  the  salinity, 
more  exactly  variation  of  the  liquid  phase  relative  content.  Such 
conditions  arise  also  at  the  ice-seawater  interface.  In  this 
connection  we  should  note  that  the  signal  which  can  be  considered 
after  reflection  from  the  lower  boundary  of  the  ice  may  begin  to 
form  somewhat  above  this  boundary,  specifically  at  the  upper  boundary 
of  the  so-called  transition  layer.  The  thickness  of  this  layer  is 
determined  by  the  temperature  conditions. 

The  electromagnetic  wave  velocity  measurements  were  made  using 
the  same  technique  and  the  results  are  shown  in  Figure  18. 

TABLE  8 


Sampling  depth,  m . . , 

. . 0-16 

17 

17.5 

18  19  22 

32 

Salinity,  % 

. . 0 
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Figure  18.  Experimental  electromagnetic  wave 
propagation  velocity  in  sea  ice  as  a function  of 
its  thickness.  1-  in  pack  ice  without  snow; 

2-  in  pack  ice  with  snow  cover  more  than  0.5  m 
thick;  3-  in  pack  ice  segment  with  porous  texture; 
4-  in  young  ice  1.5  m thick  with  air  temperature 
from  -23°C  to  25°C;  5-  ice  thickness;  6-  m/ys. 


We  did  not  make  any  special  measurements  of  the  electromagnetic 
characteristics  of  the  ice  covers  of  freshwater  bodies.  It  can  be 
stated  a priori  that  the  known  results  of  laboratory  measurements 
made  on  freshwater  ice  specimens  and  also  the  results  of  field 
measurements  made  on  glacier  ices  can  be  extended  without  major 
error  to  the  ice  covers  of  lakes,  rivers,  and  so  on.  Here,  natural- 
ly, we  must  consider  the  temperature  of  these  covers  and  introduce 
corrections  to  e ' and  tg6.  It  should  be  emphasized  that  the  physical 
characteristics  of  freshwater  body  ice  covers  are  more  stable  in 
comparison  with  the  same  characteristics  for  drifting  sea  ice.  The 
freshwater  body  ice  covers  are  more  uniform  through  their  thickness 
than  the  sea  ices,  there  are  no  salt  components  in  their  structure, 
there  is  a very  marked  ice-water  crystallization  boundary  which 
leads  to  electromagnetic  contrast  which  is  constant  and  large  in 
magnitude . 
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Specifics  and  Techniques  of  Sea  and  Freshwater  Ice  Radiosounding. 


The  radar  measurements  made  over  several  years  of  Antarctic  and 
Arctic  glacier  thickness  have  created  a good  theoretical,  procedural, 
and  engineering  basis  for  transition  to  radiosounding  of  sea  and 
freshwater  ice  covers.  Several  hardware  variants  [9]  have  been 
developed  for  radiosounding  of  floating  ice  thickness.  We  believe 
that  the  pulse  radar  method  is  promising  for  this  task.  It  is  simple 
in  technical  realization  and  from  the  viewpoint  of  reducing  the  ob- 
servational results. 

The  problem  of  remote  surveying  of  ice  thickness  includes  ob- 
taining the  above-water  and  under-water  relief  profiles.  The  survey 
scale  depends  on  the  carrier  vehicle  flight  altitude,  antenna  pattern 
width,  and  range  resolution  of  the  measurement  method  used.  In  other 
words,  the  specified  survey  scale  defines  the  radar  technical  para- 
meters and  the  aircraft  flight  conditions.  We  should  also  note  the 
influence  of  factors  such  as  scattering,  polarization,  and  focusing 
of  the  waves.  Their  evaluation  is  the  subject  of  a separate  major 
study . 

We  shall  indicate  some  general  factors  which  determine  the 
conditions  of  floating  ice  thickness  survey  with  a given  scale. 

It  is  well  known  that  the  ice  surface  and  under-water  relief 
profiles  depend  on  the  age  and  formation  conditions  of  the  ice. 
Perennial  ice  has  the  most  uneven  surface.  The  autocorrelation 
radius  L of  the  outer  surface  relief  of  perennial  ice  is  10  - 15 
meters  with  mean-square  deviation  a of  the  relief  heights  from  the 
mean  surface  equal  to  O.^to  0.5  meters.  The  autocorrelation  radius 
of  the  under-water  relief  of  perennial  ice  is  about  50  meters.  The 
maximal  ice  thickness  in  this  radius  may  reach  12-15  meters  with 
average  ice  field  thickness  of  four  meters.  For  geophysical  purposes 
it  is  sufficient  to  obtain  estimates  of  the  thicknesses  on  an  area 
with  radius  less  than  the  autocorrelation  radius  of  the  outer  surface 
relief.  In  this  case  the  thickness  estimate  will  be  biased,  depend- 
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of  the  outer  surface  relief  profile  the  following  condition  must  be 
satisfied : 

^«r.L\ 

in  which  the  left  side  is  half  the  area  of  the  first  Fresnel  zone; 
hQ  = flight  altitude;  A = radar  wavelength.  For  example,  for  radar 
carrier  frequency  6 00  MHz  resolution  of  the  outer  surface  relief  of 
ice  with  autocorrelation  radius  10  meters  can  be  obtained  with 
flight  altitude  less  than  400  meters. 


We  note  once  again  that  for  fast  remote  surveying,  the  method 
with  time  separation  of  the  signals  reflected  form  the  ice  boundaries 
is  more  convenient  from  the  technical  aspect  and  in  regard  to  inter- 
pretation of  the  results.  The  resolution  of  this  method  is  determined 
basically  by  the  radiated  pulse  duration. 


Experiments  on  radiosounding  of  the  Central  Arctic  drifting  sea 
ice  from  an  airplane  were  made  in  the  spring  of  1973  in  the  region  of 
the  SP-21  drifting  station.  A pulse  radar  with  the  following  para- 
meters was  selected  for  the  experiments : 


Carrier  frequency  100  MHz  , 

Pulse  duration 50  ns 

Antenna  directivity  6 dB 

Radar  potential 140  dB 
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The  radar  was  located  aboard  an  AN-2  airplane.  Flights  to 
measure  the  thickness  were  made  over  the  primary  ice  formations  of 
the  Central  Arctic;  perennial  pack  fields  of  thickness  more  than  3 
meters;  two-year-old  pack  of  thickness  from  2 to  3 meters;  and  young 
ice  of  thickness  1 to  1.5  meter.  The  measurements  were  made  from 
100-m  height  at  150  km/h  air  speed.  Figure  19  shows  the  pack  ice 
radar  sounding  oscillograms  taken  sequentially  in  time  with  0.5 
second  exposure.  The  time  delay  between  the  pulses  reflected  from 
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the  upper  and  lower  boundaries 
of  the  ice  are  clearly  seen  on 
the  right  side  of  the  oscillo- 
grams. The  distance  between 
their  peaks  indicates  the  pack 
ice  thickness  to  some  definite 
scale . 

A characteristic  feature 
of  the  reflected  pulses  received 
aboard  the  airplane  is  their 
short  "lifetime,"  i.e.,  the 
low  temporal  stability  of  the 
image.  The  lifetime  is  deter- 
mined by  the  carrier  airplane 
flight  speed  and  certain  other 
factors.  The  levels  of  the 
reflected  signals  recorded  in 
flight  are  determined  by  the 
state  of  the  reflecting  surfaces 
and  by  the  physical  charac- 
teristics of  the  ice.  Relief 
irregularity  and  different 
degree  of  conformance  of  the 
snow  cover  and  the  upper  ice 
layers  lead  to  fluctuations  of 
the  signal  reflected  from  the 
upper  surface  of  the  ice. 

These  same  factors  together 
with  scattering  and  absorption 
within  the  ice,  which  have 
spatial  variability,  lead  to 
fluctuations  during  flight  of 
the  signals  reflected  from  the 
lower  surface  of  the  ice. 


Figure  19.  Pack  ice  sounding 
oscillograms . 
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The  size  of  the  area  illuminated  by  the  radar  also  influences  the 
nature  of  the  fluctuations.  Signal  fluctuations  are  also  observed 
whensounding  young  ice  fields  with  smooth  boundary  surfaces,  although 
these  fluctuations  are  weaker. 

The  measurements  showed  that  the  overall  attenuation  of  the 
signals  reflected  by  the  upper  surface  of  the  ice  during  flight  above 
pack  ice  at  100  meters  altitude  varies  in  the  range  from  90  to  120  dB. 
The  theoretical  geometric  losses  are  80  dB,  the  losses  caused  by 
reflection  (Fresnel)  for  tfn  =3-5  are  10  dB,  the  losses  resulting 
from  surface  scattering  and  media  impedance  mismatch  do  not  exceed 
30  dB  (for  young  ice  the  latter  does  not  exceed  10  dB).  The  overall 
attenuation  of  the  signals  reflected  by  the  lower  surface  of  the  ice 
is  approximately  the  same  order  of  magnitude.  In  Figure  20  the 
dashed  lines  are  the  theoretical  curves  of  overall  attenuation  of 
the  signal  reflected  from  the  lower  surface  versus  ice  thickness  for 
fixed  values  of  the  specific  absorption  in  homogeneous  ice  with  plane- 
parallel  boundaries.  The  points  are  the  experimental  signal  atten- 
uation values.  We  see  from  the  figure  that  the  specific  absorption 
does  not  exceed  1-5  dB/m  for  this  series  of  studies,  conducted  above 
the  most  characteristic  Central  Arctic  ices  in  April,  1973-  The 
entire  range  of  recorded  signal  amplitudes  is  shown  in  Figure  20. 

The  signals  shown  in  Figure  19  were  recorded  with  exposure  time  0.5 
sec;  during  this  time  the  carrier  vehicle  traveled  20  meters.  There- 
fore it  seems  probable  that  statistical  interpretation  of  the  measure- 
ment results  will  be  very  promising  in  the  future.  This  requires 
photographic  recording  of  the  reflected  signals  with  short  exposure 
time  or  with  the  aid  of  a high-speed  storage  system. 

These  summertime  experiments  uncovered  still  another  important 
problem  of  ice  sounding;  namely,  the  reliability  of  the  electro- 
magnetic wave  velocities  in  the  ice.  It  is  well  known  that  the  elec- 
tromagnetic wave  velocity  in  sea  ice  is  determined  by  the  relative 
content  of  the  liquid,  solid,  and  gaseous  phases.  In  other  words, 
the  radio  signal  velocity  will  be  different  in  ices  of  different  age 
and  ice  formation  conditions.  Measurements  were  made  of  the  thickness 
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average  velocity  of  vertical 
electromagnetic  wave  propagation 
in  the  ice  fields  located  along 
the  flight  paths.  The  measure- 
ments showed  that  the  wave 
velocity  has  different  values  in 
the  same  pack  field.  On  the 
snow-free  areas  the  velocity  was 
on  the  average  150  m/usec,  on 
areas  where  the  snow  cover  ex- 
ceeds thickness  0.5  meter  the 
velocity  is  much  lower  and  de- 
creases with  reduction  of  the 
ice  thickness.  Because  of  this 
the  radar  sounding  data  will 
make  the  pack  field  appear  more  homogeneous  in  thickness  than  it 
actually  is.  In  slush  zones,  where  the  ice  has  a porous  texture, 
the  wave  velocity  is  190  - 300  m/usec,  the  wave  velocity  in  one-meter- 
thick  young  ice  with  stable  air  temperature  of  -23,  -25°C  is  90  - 
100  m/usec.  We  can  expect  that  the  rf  signal  velocity  in  young  ice 
will  also  be  determined  to'  considerable  degree  by  the  ice  temperature 
regime.  The  velocity  will  be  lower,  the  thinner  the  ice  and  the 
higher  its  temperature.  The  marked  variability  of  signal  propagation 
velocity  in  sea  ice  requires  for  correct  interpretation  of  the  radar 
thickness  survey  classification  of  the  ices  being  sounded  and  account 
for  the  thickness  of  the  snow  cover  and  the  temperatures  of  the  ice 
and  of  the  snow  cover. 

In  conclusion  we  note  that  the  series  of  studies  made  by 
Tripol'nikov  on  Central  Arctic  drifting  ice  enriched  considerably 
our  knowledge  of  the  electromagnetic  characteristics  of  the  primary 
drifting  ice  forms.  These  measurements  created  the  physical  basis 
for  radiosounding  of  sea  ice  cover  and  in  some  measure,  freshwater 
ice  cover.  The  remote  measurements  of  drifting  ice  thickness  from 
an  airplane  confirmed  the  fundamental  possibility  of  recording  sea 
ice  thickness  using  the  pulse  radar  method.  Analysis  of  these 
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Figure  20.  Overall  attenuation 
of  radio  signal  during  airplane 
sounding  of  drifting  ice  of  dif- 
ferent thickness.  1-  overall 
attenuation;  2-  ice  thickness; 

3-  dB;  dB/m. 
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experiments  disclosed  new  and  very  important  aspects  of  the  general 
problem  of  remote  measurement  of  sea  ice  cover  thickness:  the  need 

for  taking  into  account  the  influence  of  snow  cover  and  temperature 
distribution  in  the  ice  body  on  radio  signal  propagation  velocity, 
and  the  desirability  of  performing  statistical  analysis  of  the 
recorded  radio  signals  reflected  by  the  bounding  surfaces  of  the 
ice  being  measured. 

The  experiments  on  radiosounding  of  freshwater  body  ice  covers 
will  probably  differ  somewhat  from  those  discussed  above.  Specifical- 
ly, the  lower  absorption  of  rf  waves  in  the  ice  and  the  higher  stab- 
ility of  the  electromagnetic  wave  propagation  velocity  in  the  interior 
of  freshwater  ice  will  have  an  effect.  On  the  other  hand,  the  high 
(in  comparison  with  sea  ice)  electromagnetic  wave  velocity  in  fresh- 
water ice  makes  it  necessary  to  increase  the  .’adar  resolution  in  order 
to  measure  freshwater  and  sea  ices  which  are  comparable  in  thickness. 
Remote  stationary  measurements  of  lake  ice  about  two  meters  thick  were 
discussed  in  [21].*  Radiosounding  was  performed  at  three  frequencies: 
10  , 10,  and  440'MHz,  the  pulse  durations  were  15,  25  and  25  ns 
respectively.  The  experiments  showed  that  there  is  reliable  indication 
of  the  pulses  from  the  upper  and  lower  surfaces  of  ice  covers  which 
are  snow  covered  and  those  which  are  free  of  snow.  These  experiments 

| 

confirm  convincingly  the  promising  nature  of  radar  sounding  of  fresh- 
water ice.  ^ 

10.  APPLICATION  OF  THE  PULSE  RADAR  METHOD  IN  OTHER  GEOPHYSICAL 

OPERATIONS 

The  theoretical  questions  examined  above  relative  to  electro- 
magnetic wave  propagation  in  various  ices  and  their  experimental 
resolution  are  a serious  prerequisite  for  successful  solution  of  the 
problem  of  radar  sounding  of  other  media:  permafrost,  desert  sands 

to  detect  aquifer  strata,  and  also  fresh  water.  The  experience 
— 

• * Radiosounding  was  conducted  by  installations  located  on  the  ice. 


accumulated  in  sounding  ices  and  the  widespread  use  of  various  equip- 
ment for  this  purpos.e  makes  it  possible  to  examine  the  questions  of 
the  methods  and  techniques  for  sounding  these  media. 

Electrical  Properties  of  Permafrost,  Sands,  and  Fresh  Water  in 

the  rf  Band. 

At  the  present  time  permafrost  properties  have  been  studied  most 
completely  using  direct  current  [15].  The  data  available  on  perma- 
frost properties  at  the  radio  frequencies  [15]  do  not  provide  a pic- 
ture of  their  variability  as  functions  of  electromagnetic  field  fre- 
quency or  permafrost  structure  and  condition.  Because  of  this  we 
conducted  studies  of  the  electrical  properties  of  permafrost  materials 
on  specimens  obtained  from  certain  regions  of  the  Soviet  Arctic  [26]. 

Data  on  the  electrical  properties  of  sand  with  varying  water 
content  in  the  3-30  MHz  frequency  range,  obtained  at  AANIIy  are 
presented  in  [14]. 

An  Important  result  of  these  measurements  is  the  fact  that  the 
electrical  properties  of  dry  sand  at  30  MHz  ( e'  = 3,  Ni  = 0.104  dB/m) 
are  practically  identical  to  the  properties  of  fresh  water  at  -1°C 
[14].  We  would  expect  that  under  natural  conditions  the  properties 
of  dry  sand  would  be  similar  to  those  measured  under  laboratory 
conditions  while  the  properties  of  wet  sand  may  differ  considerably, 
particularly  with  the  presence  of  dissolved  salts. 

The  electrical  properties  of  fresh  water  have  been  studied  by 
several  authors  [4,  23];  however  the  data  obtained  (without  indica- 
tion of  the  chemical  composition  of  the  water)  differ  markedly. 

This  drawback  was  taken  into  consideration  in  our  measurements,  made 
in  Lakes  Gladyshevskoye  (Karelian  Isthmus)  and  Onezhskoye.  On  the 
basis  of  our  measurements  we  can  consider  that  e'  of  lake  water  in 
the  meter  band  is  about  81  and  varies  somewhat  with  temperature 
change.  In  this  same  band  the  absorptivity  is  about  1.5-2  dB/m. 


In  analyzing  the  electrical  properties  of  permafrost,  sand,  and 
fresh  water  we  can  note  the  following: 

1.  the  dielectric  permeability  of  these  media  does  not  vary 
significantly  in  the  frequency  band  above  10  MHz,  i.e.,  in  this  band 
these  media  do  not  have  significant  dispersion; 

2.  the  absorptivity  in  the  10-30  MHz  band  ranges  from  0.03  to 
2 dB/m,  sometimes  reaching  (for  permafrost)  10  dB/m.  In  all  these 
media  absorption  increases  with  increase  of  the  frequency. 

The  absence  of  dispersion  and  the  comparatively  low  absorptivity 
make  it  possible  to  select  a frequency  band  which  is  promising  for 
the  sounding  of  permafrost,  sand,  and  fresh  water.  This  band  includes 
the  10-30  MHz  frequencies.  In  certain  cases,  when  the  absorption  is 
low,  higher  frequencies  up  to  200-300  MHz  can  be  used. 

Calculation  of  Overall  rf  Signal  Attenuation  in  Layers  of 

Permafrost,  Sand,  and  Fresh  Water. 

The  overall  attenuation  calculation  presented  in  Section  5 can 
be  used  to  determine  the  signal  levels  in  vertical  sounding  of  other 
media.  The  values  of  the  Individual  attenuation  components  were 
obtained  by  Trepov.  We  shall  present  their  values: 

signal  attenuation  due  to  focusing  JV*  for  is  -5  dB 

for  sand;  -6  dB  to  -9  dB  for  permafrost,  and  -19  dB  for  water; 

the  scattering  by  the  surface  of  the  medium  is  deter- 
mined by  the  Fresnel  coefficient  of  the  air-medium  interface  and 
constitutes  for  passage  in  both  directions  0.7  dB  for  sand,  0.9  to 
2.5  dB  for  permafrost,  and  8.8  dB  for  water; 

the  losses  ,Maon,  caused  by  signal  depolarization  are  not 
known  and  are  subject  to  study;  in  many  cases  they  may  be  very  low; 
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'the  losses  Ntnpl  are  determined  by  the  Fresnel  coefficient 
of  the  interface  of  the  medium  being  sounded  and  the  underlying 
medium.  For  the  calculation  we  can  take  the  following  conditions  at 
the  lower  boundary: 

1.  "Reflection  in  sand  takes  place  from  the  surface  of  water 

concentration  change  from  0.1  (e'  = 4 for  f = 30  MHz)  to  0.2  (e'  = 10). 
Here  Kf  = 0.228,  /V^  = 12.8  dB. 

2.  For  permafrost  A\>Tp  varies  from  24.6  dB  (for  permafrost  with 
e’  = 8)  to  19  dB  (for  permafrost  with  e'  = 4).  In  the  calculation 

it  was  assumed  that  the  underlying  formation  is  wet  sand  (s'  = 10  at 
30  MHz). 

3.  When  sounding  a freshwater  basin  the  losses  (A^)  at  the 

water-bottom  interface  are  5.6  dB  for  rocky  bottom  (e'  = 8),  for 
sandy  bottom  with  water  concentration  in  the  sand  0.38  (s’  = 27.5  for 
f = 30  MHz)  = 11.6  dB).  Table  9 presents  rf  signal  absorptivity 

in  the  media  for  which  the  intrascopy  problem  is  discussed  in  the 
present  study. 

The  final  expressions  for  the  overall  losses  (in  dB)  for 
-=f2, A,-*- 0 have  the  form: 

1.  for  sand  Mr + 20  Ig  /»+  j 1; > 

2.  for  permafrost  20  lgA+ 14,5;M 

3.  for  water  Mx-*J/nlt  + 20  lg/H-3,8. 

Region  of  Application  of  Radar  Sounding  of  Sands,  Permafrost, 

and  Fresh  Water. 

On  the  basis  of  the  above  overall  attenuation  calculation, 
assuming  a radar  potential  for  example  of  160  - 170  dB , we  can  eval- 
uate the  sounded  layer  thickness  for  which  the  reflected  signal  ex- 
ceeds the  receiving  channel  noise  level.  These  data  are  presented  in 
Table  10  and  were  calculated  with  the  antennas  located  in  the  air  near 
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TABLE  9 

AVERAGED  ELECTRICAL  CHARACTERISTICS  OF  SOME  MEDIA 


Medium 

Velocity, 

m/us 

Absorptivity  at  30  MHz,  dB/m 

Cold  glaciers 

Warm  and  temperature  glaciers 

• 167 
167 

0.02-0,03 

0.06-0.08 

Permafrost 

' 110  ’ ’■  - 

1. 5-5.0  ' ' 

Freshwater  bodies 

33 

1. 0—2.0 

Sandy  soil: 

■n 

...  / * * 

• -4 

- dry 

V 1 13  _ 

1. 5-2.0 

- wet 

.55  ^ 

Pack  Sea  Ice 

1 150  :,i 

l— 2n  ’ 

Young  Sea  Ice 

100  1 

i 3--o  ' 

the  upper  surface  of  the  medium  (except  ->r  water). 

Technique  and  Procedural  Characteristics  of  Radar  Sounding  of 
Strongly  Absorbing;  Media. 


If  we  consider  the  primary  objective  of  radar  sounding  of  media 
to  be  the  measurement  of  the  thickness  of  their  layer  to  the  bound- 
aries with  the  underlying  formations  or  to  the  temperature  jump 
surfaces,  then  by  analogy  with  glacier  sounding  we  can  identify  three 
types  of  problems,  the  technique  and  procedure  for  the  solution  of 
which  differ: 

1.  measurement  of  the  maximal  layer  thicknesses; 

2.  measurement  of  thicknesses  which  are  definitely  less  than 
the  maximal  thicknesses  measurable  with  the  aid  of  radar  sounding; 

3.  measurement  of  small  thicknesses,  for  example  the  freezing 
depth  of  soil  or  of  thin  layers  of  sand  and  water. 
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THICKNESS  OP  SOUNDED  LAYER  FOR  SOME  MEDIA 
(RADAR  POTENTIAL  160  - 170  dB) 


MEDIUM 

RADAR  CARRIER 
FREQUENCY,  MHz 

ICE  THICKNESS,  M 

Dry  Sand 

30 

600 

Wet  Sand: 

- Water  concentration  0.1  to 

30 

140 

Without  salts  0.2  to 

30 

40 

Rock-ice  permafrost  (t=-15°C) 

10 

250 

Frozen  soil 

10 

15 

Fresh  water: 

- from  [17] 

30 

100 

- from  [4] 

30 

3000 

The  selection  of  equipment  parameters  such  as  carrier  frequency, 
radiated  power,  frequency  spectrum  width  (pulse  duration  for  pulse 
radars),  and  processing  in  the  receiver  do  not  present  any  difficulty. 
Typical  oarameters  of  the  radar  for  measuring  the  maximal  thickness 
of  permafrost,  sand,  and  fresh  water  may  be  as  follows: 


Carrier  frequency 
Pulse  power 
Pulse  duration 
Receiver  bandwidth 
Radar  potential 

The  parameters  of  the  radar 
type  may  be  as  follows: 

Carrier  frequency 
Pulse  power 
Pulse  duration 
Receiver  bandwidth 
Radar  potential 


10  MHz 
100  kW 
0.3  us 
3 MHz 
183  dB 

for  solving  problems  of  the  second 


30  MHz 
10  kW 
0.1  us 
15  MHz 
165  dB 


76 


.IV-** 


1 


i 

* 

i 

V 

f 

( 

* 

5 


f 

« 


5 

I 


Typical  parameters  of  the  radar  for  measuring  thin  layers  of 
permafrost,  sand,  and  water  are: 

Carrier  frequency 

Pulse  power 

Pulse  duration 

Receiver  bandwidth 

Antennas 

Radar  potential 

In  spite  of  the  fact  that  radar  sounding  of  the  considered  media 
is  possible  and  suitable  apparatus  can  be  developed,  in  the  practical 
application  of  the  method  there  are  several  serious  difficulties: 

1.  The  high  noise  level  in  the  form  of  reflections  from  irreg- 
ularities on  the  surface  of  the  medium  and  nonhomogeneities  within 
the  medium  may  make  it  difficult  to  identify  and  interpret  the  useful 
signal ; 

2.  Noise  in  the  form  of  direct  signal  transmission  from  the 
transmitter  output  to  the  receiver  input  and  reflections  from  the 
surface  increase  the  minimal  measurable  thickness.  Considering  that 
the  absolute  delay  time  for  the  maximal  measurable  thickness  is 
short  (2  to  10  us),  the  thickness  range  which  can  be  measured  in 
practice  may  be  too  small. 

We  shall  not  discuss  the  special  procedural  and  hardware  tech- 
niqhes  which  can  be  used  to  eliminate  these  problems.  We  shall  show 
that  two  of  the  problems  posed  above,  namely  radiosounding  of  fresh 
water  and  sand,  have  been  resolved  experimentally. 

First  Experiments  on  Radar  Sounding  of  Fresh  Water  and  Sand. 

The  first  radar  sounding  of  fresh  water  was  accomplished  in  Lake 
Gladyshevskoe  in  August  1969-  A pulse  radar  fabricated  by  Radiophysics 
Division  of  AANII  was  used  for  sounding.  Its  specifications  were: 
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150  MHz 
100  W 
0.015  us 
70  MHz 
wideband 
120  dB 
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Carrier  frequency 

Pulse  power 

Pulse  duration 

Pulse  repetition  frequency 

Receiver  passband 

Potential 


45  MHz 


0.1  us 
100  Hz 
10  MHz 
110  dB 


The  receiving  and  transmitting  antennas  — halfwave  dipoles  with 
reflectors  in  the  form  of  metal  sheets  — were  submerged  in  the  water. 
In  the  sounding  process  we  recorded  the  signals  from  an  artificial 
reflector  (metal  sheet)  and  from  the  bottom  of  the  lake  for  various 
depths  up  to  the  maximal  depth  for  this  particular  lake  of  15  meters. 
The  lake  botton  consisted  of  sand,  mud,  and  rocky  strata. 

The  experiments  with  the  artificial  reflecter  made  it  possible 
to  determine  the  radio  pulse  velocity  in  the  water;  it  was  33-3  m/us 
(e1  = 8l).  We  also  obtained  from  the  sounding  results  the  specific 
attenuation  of  the  radio  signal;  for  this  particular  body  of  water 
.Vi,  was  equal  to  1.4'  dB/m. 

In  1971  studies  were  made  in  the  Petrozavodsk  Inlet  of  Lake 
Onezhskoye.  In  these  experiments  we  used  a more  powerful  radar  with 
carrier  frequency  60  MHz  and  pulse  duration  0.4  usee.  The  radio- 
sounding  was  conducted  with  ship  speed  about  10  km/h.  The  antennas 
were  submerged  in  the  water.  The  radio  pulse  propagation  velocity 
in  this  lake  was  33  m/usec  ( e ’ = 82.6)  with  water  temperature  10°C. 

The  specific  absorption  was  different  on  the  different  tacks:  2 and 

1.75  dB/m.  High  values  of  the  specific  resistance  (2  dB/m)  were 
obtained  on  tacks  passing  near  urban  sewer  outfalls.  It  was  noted 
that  the  different  geological  materials  forming  the  bottom  have 
different  reflectivity.  The  values  of  tg5  were  calculated  from  the 
experimental  data  for  the  waters  of  Lakes  Onezhskoye  and  Gladyshevskoye 
and  compared  with  the  results  of  laboratory  measurements  performed 
by  Saxton  and  other  investigators.  Radiosounding  of  Lakes  Glady- 
shevskoye and  Onezhskoye  shows  that  the  radar  method  makes  it  pos- 
sible to  measure  comparatively  simple  the  depth  of  fresh  water  bodies, 
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detect  objects  in  the  water,  determine  the  electrical  properties 
of  the  water  under  natural  conditions,  and  evaluate  the  degree  of 
contamination  of  water  bodies. 

In  July  1972  radiosounding  of  gravelly  sand  was  carried  out  by 
the  Valdai  Scientific  Research  Hydrological  Laboratory  of  the  State 
Hydrological  Institute  in  order  to  determine  groundwater  depth. 
According  to  hydrological  observation  data  groundwater  lies  at 
depths  of  6.7  and  13  meters  in  the  regions  studied.  Sounding  was 
performed  with  a radar  with  the  following  characteristics: 


Carrier  frequency 
Pulse  duration 
Radiated  power 
Potential 


60  MHz 
0.1  ys 


150  dB 


Several  series  of  radiosoundings  were  conducted,  differing  in 
mutual  positioning  of  the  receiving  and  transmitting  antennas  (half- 
wave dipoles),  field  polarization,  and  distance  between  the  antennas. 
To  facilitate  interpretation  of  the  received  signals  the  antennas 
were  positioned  in  the  upper  layer  of  the  sand  at  depths  from  0.1  to 
0.4  meters.  We  also  used  reflectors  in  the  form  of  sheets  and  tubes 
positioned  on  the  sand  surface  above  the  antennas.  The  experiments 
disclosed  the  presence  of  four  characteristic  signals:  direct  signal 

from  antenna  to  antenna,  traveling  through  the  air;  direct  signal  from 
antenna  to  antenna,  traveling  in  the  ground;  signal  scattered  by  the 
upper  surface  of  the  water-saturated  layer;  and  in  certain  cases  the 
signal  scattered  by  the  lower  surface  of  this  layer.  Analysis  of  the 
oscillograms  showed  that  in  the  region  with  the  aquifer  layer  located 
at  a depth  of  6.7  m the  signal  reflected  from  the  layer  lying  at  the 
7 m depth  is  noted:  the  layer  thickness  was  about  3 m.  In  the  second 

case,  when  the  layer  was  located  at  a depth  of  13  m,  radiosounding 
also  showed  this  depth,  i.e.,  13  m.  The  thickness  of  the  water- 
saturated  layer  from  the  radiosounding  data  was  11  m in  this  case. 


This  sounding  investigation  is  believed  to  be  the  first 
experimental  proof  of  the  feasibility  of  searching  for  and  locating 
groundwater  with  the  use  of  radar  techniques. 
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